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This  report  represents  the  results  of  the  efforts  expended  in  performance 
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Bell  Helicopter  Company  Project  Engineer.  Mr.  Robert  Nicholson  was 
the  Air  Force  Project  Engineer. 

This  final  report  is  presented  in  four  volumes.  The  first  describes  the 
mathematical  model  and  the  methods  used  to  calculate  stability  character¬ 
istics.  They  are  of  sufficient  complexity  that  a  digital  computer  is 
necessary  for  the  solution  of  the  equations.  The  second  volume  presents 
the  results  of  sample  computations  and  discusses  input  and  output  formats 
and  potential  problems.  In  the  third  volume,  the  computer  program  cross- 
reference  indices  are  given  with  explanation  of  FORTRAN  variable  and 
subroutine  names.  Volume  IV  is  a  complete  FORTRAN  listing  of  the  computer 
program. 
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ABSTRACT 


This  report  describes  a  mathematical  model  of  rotorcraft  that 
may  be  used  to  determine  characteristics  of  performance,  sta¬ 
bility,  response,  and  rotor  blade  loads.  The  complexity  of 
the  equations  used  requires  the  use  of  a  digital  computer  for 
efficient  solution.  This  four-volume  report  describes  the  com¬ 
puter  program  in  detail  and  illustrates  the  method  of  computing 
rotorcraft  characteristics  by  specific  example. 

This  volume  presents  an  overview  of  the  computer  program  capa¬ 
bilities  and  describes  the  formulation  of  the  mathematical 
model.  The  second  and  third  volumes  present  sample  cases, 
detailed  input  and  output  formats,  and  cross-reference  indices 
of  FORTRAN  variable  and  subroutine  names  and  usage.  A  complete 
FORTRAN  listing  comprises  all  of  the  fourth  volume. 

The  computer  program  was  modified  from  an  earlier  version 
so  that  the  stop-fold-rotor  concepts  could  be  represented.  The 
two  concepts  considered  are  the  tilt-forward-trail-aft  (TFTA) 
and  the  horizontal -stop-fold  (HSF)  configurations.  The  TFTA 
aircraft  simulation  uses  wing-tip  mounted  rotors  for  helicopter 
flight  which  rotate  through  90  degrees  to  act  as  propellers 
in  airplane  flight.  If  auxiliary  propulsion  is  available, 
these  rotors  may  be  feathered,  stopped,  and  folded  back  for 
high-speed  flight.  The  HSF  helicopter  simulation  uses  a  con¬ 
ventional  main  and  tail  rotor  arrangement.  The  addition  of  a 
high-capacity  rotor  brake  enables  the  rotors  to  be  stopped  in 
flight  at  a  specified  location  and  then  folded  back  through  a 
predetermined  schedule.  For  both  of  these  aircraft,  stability 
and  control,  performance,  and  response  characteristics  may  be 
computed  at  any  time  during  the  conversion  maneuvers . 

No  small  angle  assumptions  are  made  in  the  analysis.  Aero¬ 
dynamic  interference  effects  between  rotors,  wings,  elevators, 
and  fins  have  been  represented.  First-order  elastic  effects 
of  the  fuselage  and  wing  have  been  included  also.  Partial  and 
total  stability  derivatives  may  be  computed  in  accelerated  or 
steady  flight  conditions.  Root  locations  of  the  rigid-body 
aircraft  motion  are  computed  as  well  as  the  transfer  functions 
of  the  pilot  controls.  The  effects  of  automatic  controls, 
gusts,  and  weapon  recoil  may  be  easily  computed  in  the  form  of 
time  histories  of  motion. 
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SECTION  I 


INTRODUCTION 

Because  the  equations  that  describe  Lhe  physical  phenomena  are 
quite  complex,  and  because  the  mathematical  model  is  designed 
to  represent  a  wide  variety  of  rotorcraft  configurations,  stabi 
lity  and  control  characteristics  are  calculated  by  a  digital 
computer.  The  mathematical  model  has  evolved  from  a  relatively 
simple  rotor-loads  analysis  over  some  ten  years.  As  it  was 
refined,  the  program  was  carefully  checked  at  each  step  to  en¬ 
sure  that  it  included  all  the  modifications  of  the  mathematical 
model.  The  original  rotor-loads  analysis  emphasized  dynamics; 
later,  aeroelastic  characteristics  were  added.  The  capabili¬ 
ties  of  the  analysis  have  been  demonstrated  in  published  papers 
with  some  sample  cases  that  show  the  correlation  with  measured 
data. (1*2)  Adding  the  aerodynamics  of  the  wing,  fuselage,  and 
empennage  to  the  mathematical  model  of  the  rotor  made  the 
representation  of  trimmed  flight  conditions  easier  and  more 
accurate.  Then  another  section  was  added  to  the  program  to 
enable  it  to  represent  maneuvers.  It  then  became  possible  to 
compute  the  blade  loads  in  maneuvers  as  well  as  in  trimmed 
flight . 

In  1965,  Bell  entered  into  a  contract  with  the  U.S.  Army 
Aviation  Materiel  Laboratories  to  determine  analytically  the 
gust  response  of  rotors.  (3)  The  maneuver  section  of  the 
analysis  was  augmented  during  the  performance  of  this  contract 
to  compute  response  to  gusts  of  any  magnitude,  from  any  direc¬ 
tion.  The  gust  shapes  could  be  sharp-edged,  ramp,  or  sine- 
squared.  A  version  of  the  computer  program  was  delivered  to 
the  Army  under  this  contract,^'  and  a  subsequent  version^)  is 
now  operational  on  the  CDC6600  digital  computer  at  the  NASA 
Langley  Research  Center. 

The  analysis  was  further  expanded  during  studies  of  composite 
research  aircraft  to  include  tilting  proprotors  and  a  more 
sophisticated  control  system. (6)  Concurrently,  a  subroutine 
that  enabled  the  computation  of  partial  stability  derivatives 
in  trimmed  flight  was  added.  Subsequent  modifications  of  the 
trim  analysis  included  the  equations  which  represent  accelera¬ 
ted  flight  conditions:  banked  turns,  pull-ups,  or  push-overs; 
during  level,  climbing,  descending,  accelerating,  or  decelera¬ 
ting  flight.  The  program  can  now  represent  a  variety  of  air¬ 
craft  configurations,  including: 

Helicopters 

Single  Main  and  Tail  Rotors 

Tandem  Rotors 

Laterally-Displaced  Rotors 

Compound  (jet  or  propeller  propulsion) 
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-  Fixed-Wing  Aircraft 

Jet-Propelled 

Propeller-Driven,  Single  or  Twin  Engines 

-  Composite  Aircraft 

Tilt -Rotor 

Tilt-Stop -Fold -Rotor 
Edgewise-Stopped-Rotor 

The  maneuver  analysis  can  also  be  used  to  evaluate  the  effects 
of  a  wide  variety  of  inputs  and  systems  on  aircraft  response. 
Examples  of  the  kinds  of  inputs  it  can  handle  are  gusts,  weapon 
recoil,  drop  of  external  stores,  control  motions,  stability 
augmentation  signals,  and  the  commands  of  rotor-pitch  governors. 

Sample  cases  of  stop-fold  maneuvers  have  been  computed  and  the 
resulting  time  histories  are  discussed  in  Volume  II,  Section 
IID.  These  cases  are  illustrative  only  and  do  not  represent 
any  particular  aircraft  configuration. 
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SECTION  II 


COMPUTER  PROGRAM  ORGANIZATION 


The  computer  program  is  divided  into  several  major  sections. 
The  relationships  between  the  major  sections  are  illustrated 
by  the  flow  chart  in  Figure  1.  The  sections  of  the  program 
which  are  used  in  any  particular  problem  are  controlled  by 
the  input  data  cards  for  that  section.  Therefore,  the  number 
of  cards  of  input  data  required  can  vary  considerably,  depend¬ 
ing  upon  how  the  program  is  to  be  used.  After  the  data  cards 
are  read,  problem  constants  are  computed  and  parameters  are 
initialized.  The  functions  of  the  subsequent  computational 
blocks  of  Figure  1  are  discussed  below. 

A.  DETERMINATION  OF  EQUILIBRIUM  FLIGHT  CONDITION 

The  six  conventional  equations  describing  aircraft  force  and 
moment  equilibrium  must  be  satisfied  for  the  trim  flight  con¬ 
dition.  Additionally,  the  equations  describing  the  longitu¬ 
dinal  and  lateral  moments  acting  on  both  rotors  must  be 
satisfied.  Since  the  rotor  speed  is  assumed  constant,  there 
are  ten  equations  to  be  satisfied. 

The  ten  unknowns  to  be  determined  are:  the  longitudinal  and 
lateral  flapping  angles  of  each  rotor,  the  pilot  contro1 
positions,  and  the  angular  orientation  of  the  aircraft  In 
space.  The  pilot  controls  are  the  collective  pitch  lever, 
the  longitudinal  and  lateral  cyclic  stick,  and  the  pedals. 

The  Euler  angle  of  pitch  is  always  determined  by  the  trim 
technique  and  either  the  roll  or  yaw  angle  is  computed  as 
determined  by  the  input  data.  The  mathematical  techniques 
used  in  the  trim  subroutine  are  discussed  in  detail  in 
Section  IV.  Essentially,  the  trim  problem  is  one  of  deter¬ 
mining  the  solution  of  ten  nonlinear  equations  in  ten  un¬ 
knowns  . 

The  equilibrium  flight  condition  may  be  any  one  of  three 
conditions  described  below.  For  all  conditions,  the  rate  of 
climb  or  descent  must  be  specified. 


1.  One-g,  Wings-Level  Flight 

The  first,  and  most  common,  condition  is  that  of  unaccelerated 
wings-level  flight.  Trim  parameters  are  computed  at  a  speci¬ 
fied  roll  or  yaw  angle.  At  speeds  below  about  50  knots  for 
helicopters,  the  yaw  angle  must  be  specified  (usually  zero)  so 
equilibrium  is  determined  by  means  of  roll  angle  orientation 
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of  the  gravity  vector.  Generally,  only  a  few  degrees  of  roll 
are  required.^  If  the  roll  angle  is  set  to  zero  at  low  speeds, 
the  aerodynamic  forces  in  the  yaw  moment  equation  are  so  small 
that  very  large  yaw  angles  may  be  required  to  trim.  It  is 
possible  that  no  trim  solution  can  be  found  with  zero  roll 
angle  input  at  low  speeds.  At  speeds  higher  than  about  50 
knots,  either  roll  or  yaw  may  be  specified. 

The  determination  of  trim  conditions  at  speeds  over  150  knots 
can  be  difficult  because  of  the  increasingly  nonlinear  nature 
of  the  trim  equations  caused  by  rotor  stall,  compressibility, 
and  reverse  flow.  The  only  effective  method  for  overcoming 
the  high  speed  trim  problem  to  date  has  been  the  use  of  the 
parameter  sweep  feature  of  the  program  to  obtain  trim  condi¬ 
tions  at  successively  higher  speeds.  The  starting  trim  para¬ 
meters  at  each  speed  are  the  final  trim  values  at  the  pre¬ 
ceding  lower  speed.  By  reducing  the  change  in  speed  from 
20  knots  to  10,  5,  or  even  fewer  knots,  either  a  higher  speed 
trim  point  is  obtained,  or  the  reasons  trim  cannot  be  obtained 
will  usually  be  apparent  after  examination  of  the  preceding 
trim  parameters  and  their  variation  with  speed. 

Probably  the  most  common  source  of  trouble  is  that  the  power 
increases  beyond  the  capability  of  the  anti-torque  system  of 
conventional  helicopters  to  counteract.  In  this  case,  the 
pedal  control  exceeds  the  travel  limits  and  directional  sta¬ 
bility  is  reduced  because  of  tail  rotor  stall.  The  user  must 
then  do  what  the  pilot  does  and  begin  (simulate)  a  descent  on 
the  computer.  Simple  energy  theory  provides  a  fairly  accurate 
estimate  of  the  rate  of  descent  required. 

R/D  =  5|0  (HPrequired  -  HPavaUable)  ft/sec 
2 .  Push-overs  and  Pull-ups 

The  trim  parameters  may  be  determined  during  symmetrical 
(wings-level)  push-overs  and  pull-ups.  Since  this  maneuver 
is  not  one  that  can  be  maintained  indefinitely,  the  computed 
data  are  just  a  ’'snapshot*'  of  the  flight  conditions  existing 
at  a  specified  instant.  Rate  of  descent  in  the  actual  maneuver 
is  never  constant,  for  example,  but  trim  conditions  are  com¬ 
puted  at  the  rate  of  descent  specified. 

The  accurate  simulation  of  this  maneuver  in  the  trim  subroutine 
is  complicated  by  the  capability  of  the  pilot  of  uie  helicopter 
to  induce  normal  accelerations  by  txv*  controls:  the  collective 
pitch  lever  and  the  longitudinal  cyclic  stick.  The  normal  trim 
technique  used  in  the  computer  is  to  adjust  these  controls 
simultaneouslv  ^  attain  V  e  desired  flight  condition.  This 
may  not  be  tvy  *ay  the  pilot  actually  controls  the  aircraft, 
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depending . upon  speed  and  the  amount  and  direction  of  normal 
acceleration  desired.  Push-over  maneuvers  induce  less  than 
one-g  loading  and  are  usually  performed  by  the  pilot  with  the 
collective  pitch  lever  fixed.  Pull-up  maneuvers  are  usually 
done  with  the  collective  pitch  lever  fixed  if  the  higher  g 
loading  will  not  be  maintained  very  long,  that  is,  for  small 
flight  path  corrections,  or  if  altitude  loss  is  not  signifi¬ 
cant.  If  these  conditions  are  not  true,  then  the  pilot  tends 
to  use  collective  and  cyclic  control  together.  Both  of  these 
control  techniques  may  be  simulated  in  the  trim  subsection. 
Collective  pitch  is  used  alone  at  very  low  speeds  and  in  hover 
to  change  g  loading.  This  is  most  easily  simulated  by  varying 
the  gross  weight  for  which  trim  parameters  are  being  computed. 
Division  by  the  actual  gross  weight  then  gives  the  desired 
g-loading. 

To  simulate-  simultaneous  control  usage,  the  roll  angle  is  set 
to  zero  and  the  g-load  required  is  specified  in  the  input  data. 
The  normal  trim  techniques  are  used.  Trim  conditions  at  less 
than  one-g  have  always  been  determined  easily  by  this  techni¬ 
que.  At  higher  loadings  than  about  1.5  g,  the  above  trim 
technique  will  fail  to  converge  to  a  solution.  The  reasons 
for  this  failure  to  converge  have  not  been  determined  to  date 
but  initial  work  indicates  that  this  occurs  when  the  main 
rotor  power  exceeds  the  available  power.  In  other  words,  one 
can  not  simulate  on  the  computer  a  flight  condition  which  is 
not  attainable  in  flight  and  expect  a  solution  to  be  found. 
There  are  two  ways  to  obtain  the  solution  at  higher  g  levels; 
however,  both  involve  a  loss  of  energy.  If  airspeed  must  be 
maintained  (kinetic  energy),  then  a  rate  of  descent  must  be 
used  (loss  of  potential  energy)  to  reduce  the  main  rotor  power 
requirement  and  enable  higher  normal  accelerations  to  be 
reached.  The  other  method  involves  loss  of  airspeed  and  no 
change  of  altitude. 

Since  the  trim  equations  assume  constant  airspeed,  it  is 
necessary  to  simulate  a  longitudinal  acceleration,  that  is, 
a  force,  acting  through  the  center  of  gravity  which  is  equiva¬ 
lent  to  a  power  through  the  familiar  energy  equation. 

Force  -  vei0City  ^HPrequired  "  HPavailable^ 

This  is  done  by  altering  the  input  data  in  the  following 
manner: 

-  Main  rotor  collective  pitch  is  input  directly  and 
the  rotor  collective  pitch  is  locked.  The  value  used 
may  be  the  one-g  value,  or  the  highes t-attainable-g 
val ue . 


5 


A  single  auxiliary  jet  is  located  at  the  center  of 
gravity  and  oriented  along  the  approximate  flight  path 
at  the  highest -attainable -g  condition.  Little  error 
results  if  it  is  simply  oriented  along  the  X-axis. 

-  The  range  of  jet  thrust  is  specified  so  that  all 
horizontal  forces  likely  to  be  needed  are  available. 

The  pilot's  collective  pitch  lever  is  linked  to  the 
auxiliary  jet  control. 

Once  this  is  done  the  usual  set  of  trim  equations  is  solved, 
and  the  required  value  of  jet  thrust  is  determined.  This 
force,  when  divided  by  the  weight,  is  representative  of  the 
linear  deceleration  that  would  be  acting  on  the  aircraft  in 
such  a  high-g  flight  condition. 

3 .  Banked  Turns 


Trim  parameters  may  be  computed  for  turning  flight.  The  normal 
acceleration,  bank  angle,  or  turn  radius  may  be  specified  in 
the  input  data.  Assuming  that  the  velocity  is  known,  g  is  the 
acceleration  of  gravity,  and  either  n,  <J> ,  or  R  are  specified; 
the  unspecified  quantities  are  determined  by  matrix  transforma¬ 
tion  of  the  Euler  angles  and  velocities.  For  level  turns  the 
transformation  matrix  may  be  reduced  to  the  following  equations, 
assuming  that  yaw  angle  is  zero. 


n  =  sec<|> 


The  angular  rates  of  the  aircraft  about  the  X,  Y,  and  Z  axes 
are : 

p  =  -  &  sine  Ai2  -  1 
q  =  ^  cosq  (n2  -  1) 
r  =  — ^  cosq  >/n2  -  1 

These  angular  rates  and  normal  accelerations  are  used  in  the 
aircraft  equations  of  motion,  which  are  discussed  in  Section 
IIIG,  to  determine  the  inertial  forces  required  to  accurately 
simulate  the  flight  conditions. 


Coordinated  banked  turns  always  induce  more  than  one-g  load¬ 
ings,  of  course.  The  same  difficulties  in  obtaining  the  trim 
solution  to  high-g  turns  are  encountered  as  previously  dis¬ 
cussed  for  symmetrical  pull-ups.  The  same  means  of  obtaining 
solutions  may  be  used  as  described  previously. 
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B.  STABILITY  ANALYSIS 

After  determination  of  the  trim  parameters,  the  stability 
derivatives  may  be  evaluated  about  the  trim  point  if  desired. 
These  derivatives  are  the  total  aircraft  derivatives  with 
respect  to  linear  and  angular  velocities  and  are  evaluated 
in  the  body  axis  system.  The  technique  used  is  to  change  the 
trim  flight  values  of  linear  and  angular  rates  one  at  a  time, 
by  a  small  amount  which  is  specified  by  the  input  data.  After 
a  trim  value  is  incremented,  the  rotor  is  allowed  to  re¬ 
establish  flapping  equilibrium  with  the  attendant  changes  in 
forces  and  moments.  The  forces  and  moments  in  other  aero¬ 
dynamic  surfaces  also  change.  The  trim  flight  values  of 
forces  and  moments  r re  then  subtracted  from  the  incremented 
values,  and  the  differences  are  then  divided  by  the  amount 
of  the  increment  used.  The  result  is  the  total  derivative 
of  force  ot  moment  with  linear  or  angular  velocity. 

This  technique  of  evaluating  stability  derivatives  is  accurate 
as  long  as  the  increment  size  is  properly  chosen.  If  too 
small  an  increment  is  used,  then  errors  are  caused  by  small 
differences  of  large  numbers.  If  too  large  an  increment  is 
used,  then  errors  are  caused  by  the  non-linear  character  of 
many  of  the  parameters.  That  is,  the  slope  of  the  line  be¬ 
tween  two  points  which  are  too  far  apart  differs  significantly 
from  the  tangent  to  the  curve  at  the  trim  point.  The  effect 
of  increment  size  on  several  primary  derivatives  of  longitu¬ 
dinal  motion  is  illustrated  in  Figure  2.  The  values  of  the 
derivatives  are  indicated  as  are  the  resulting  root  locations 
of  the  short  period  and  phugoid  mode.  The  STAB  increment  size 
is  the  value  used  in  the  input  data.  The  correct  linear  velo¬ 
city  increment  is  about  2%  of  the  forward  velocity.  The  magni¬ 
tude  of  the  angular  rate  increment  is  always  one-hundredth 
the  linear  rate  increment,  and  its  dimensions  are  radians  per 
second.  For  example,  if  the  forward  velocity  was  200  ft/sec, 
the  linear  rate  increment  should  be  4  ft/sec  and  the  angular 
rate  increment  would  then  be  .04  rad/sec. 

The  detailed  discussion  of  the  output  format  of  the  stability 
analysis  is  given  in  Volume  II,  Section  IIIH.  A  sample  of  this 
output  is  given  in  Table  I.  The  first  force  and  moment  array 
is  the  trim  solution.  These  numbers  represent  the  force 
(pounds)  and  moment  (pound-feet)  contributions  of  the  indivi¬ 
dual  components  of  the  aircraft  to  the  body-axes  forces  and 
moments.  The  second  array  shows  the  forces  and  moments  after 
the  w-velocity  has  been  incremented.  The  third  array  is  the 
difference  between  the  first  two  arrays  and  illustrates  the 
contribution  of  each  component  to  the  w  derivative  of  interest. 
These  data  are  printed  out  for  the  three  linear  velocities, 
u,  v,  and  w,  as  well  as  the  three  angular  velocities,  p,  q, 
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and  r.  All  of  the  resulting  derivatives  are  summarized 
for  the  whole  aircraft  and  for  each  rotor.  An  example 
of  these  data  is  given  in  Table  II.  The  units  of  the 
numbers  of  Table  II  are  pounds,  pound-feet,  feet  per  sec¬ 
ond,  and  radians  per  second,  as  appropriate. 

Not  all  of  the  computed  stability  derivatives  are  currently 
used  in  the  present  analysis.  The  present  analysis  uses  the 
standard  uncoupled  equations  of  motion  which  are  available  in 
textbooks . i *2)  The  appropriate  derivatives  of  Table  II  are 
substituted  and  the  roots  of  the  characteristic  equation  are 
determined  as  are  the  phase  and  amplitude  relationships  defin¬ 
ing  the  mode  shapes. (i2)  The  output  data  computed  for  the 
longitudinal  and  lateral-directional  modes  of  motion  are 
illustrated  in  Table  III. 

The  derivatives  of  control  motion  with  respect  to  forces  and 
moments  are  obtained  from  the  partial  derivative  matrix  of 
ten  equations  in  ten  unknowns  which  is  used  to  obtain  the  trim 
solution.  These  derivatives  are  also  used  to  obtain  the 
numerators  of  the  transfer  functions  indicated  in  Table  III. 
These  data  are  often  useful  in  the  determination  of  stability 
augmentation  system  parameters. 

C.  MANEUVERS 

The  equations  of  motion  listed  in  Table  X  and  discussed  in  Sec¬ 
tion  II IG  may  be  solved  in  the  maneuver  section  of  the  program 
while  control  motions,  automatic  systems,  or  external  distur¬ 
bances  may  be  introduced  and  their  effects  computed.  This 
section  of  the  computer  program  uses  the  four-cycle  Runge-Kutta 

technique^14^  for  determining  the  values  of  parameters  at 
the  next  point  in  time.  This  section  of  the  program  provides 
full-coupled,  nonlinear  response  data  in  contrast  to  the  un¬ 
coupled,  linear  data  computed  in  the  stability  section.  The 
Runge-Kutta  technique  is  described  briefly  as  follows: 

1.  Accelerations  at  t  =  are  computed  based  upon 
the  positions  and  velocities  at  t^ . 

2.  Using  the  values  from  step  1,  the  positions  and 
velocities  are  computed  at  t  =  t^  +  1/2  At. 

3.  Using  the  positions  and  velocities  of  step  2,  the 
accelerations  at  t  =  t^  +  1/2  At  are  computed. 
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4.  Using  the  positions  and  velocities  of  step  1»  and  the 
accelerations  of  step  3,  the  new  positions  and  velo¬ 
cities  at  t  =  t^  +  1/2  At  are  computed . 

5.  Using  the  positions  and  velocities  of  step  4, 

another  acceleration  at  t  =  t,  +  1/2  At  is 
computed .  1 

6.  Using  the  positions,  velocities,  and  accelerations 
of  steps  4  and  5,  the  positions  and  velocities  at 
t  =  t^  +  At  are  computed . 

7.  Using  the  positions  and  velocities  of  step  6,  the 
accelerations  at  t  =  t^  +  At  are  comr  ted. 

8.  The  final  positions  and  velocities  t  =  t^  +  At 

are  given  by  the  following  equatio  ,  where  P  re¬ 
presents  position,  V  represents  ve.jcity,  A  repre¬ 
sents  acceleration,  and  the  subscript  denotes  which 
of  the  above  steps  is  used  to  obtain  the  value. 

Pt+At  =  P1  +  £  [V1  ♦  V6  +  2  <V2  +  V] 

Vt+At  =  VL  *  ¥  [A1  +  A7  +  2  (A3  +  A5>] 

1 .  Time  Step  Size 

The  mathematical  model  of  the  rotor  changes  significantly  in 
the  maneuver  section.  The  trim  and  stability  sections  use 
twelve  azimuth  locations  to  compute  airloads.  The  maneuver 
section  uses  as  many  azimuth  stations  as  there  are  rotor  blades. 
This  fact  largely  determines  the  time  step  that  the  user  must 
specify  in  the  input  data.  This  time  step  may  be  changed  dur¬ 
ing  a  maneuver  in  order  to  reduce  computing  time.  The  time 
step  is  determined  by  the  following  consideration  which  relates 
to  the  accurate  solution  of  the  rotor  blade  equation  of  motion. 

The  time  step  *•  ist  be  chosen  so 
that  neither  rotors  rotate  through 
more  than  about  45  degrees  of 
azimuth  during  the  time  step. 

Since  conventional  helicopters  have  tail  rotors  which  turn 
at  about  five  times  the  speed  of  the  main  rotor,  the  time 
step  is  usually  based  upon  tail  rotor  speed.  For  example, 
the  UH-1D  uses  a  main  rotor  speed  of  324  RPM  at  which  the 
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tail  rotor  speed  is  1643  RPM.  Since  45°  is  one-eighth  of  a 
revolution  then  the  largest  allowable  time  step  is  approxi¬ 
mately  60/1643(8)  =  .00455  seconds;  which  could  probably  be 
rounded  off  to  .005  seconds  without  encountering  difficulty. 

If  the  maneuver  being  performed  involved  slowing  the  rotor, 
then  the  time  step  could  be  increased  proportionately  as  the 
rotor  speed  decreased. 

2 .  Output 

The  output  format  is  detailed  in  Volume  II,  Section  III. 

Here  will  be  described  briefly  what  type  of  data  are  available 
to  the  user  from  this  section.  There  are  two  basic  types  of 
output  data:  the  maneuver  data  page  which  is  printed  at  every 
nth  time  point,  where  n  is  specified  in  the  input  data,  and 
the  plots  of  variables  versus  time  for  every  mth  time  point, 
where  m  is  specified  in  the  input  data. 

There  are  up  to  266  variables  which  may  be  plotted  versus 
time.  These  variables  relate  to  nearly  all  combinations  of 
the  following: 

(linear  1  (position  (rotor  tip  path  planes 

angular/  <velocity  1  ~  rotor  blades  (individual) 

[acceleration  all  controls 

l fuselage 

fuselage 
wing 
fin 

elevator 
jets 

weapon  recoil 

Appropriate  scaling  factors  for  the  variables  chosen  are 
determined  by  the  user  and  up  to  three  variables  may  be 
plotted  in  the  same  graph.  The  plotting  routine  used  is  a 
digital  one  which  uses  ten  inches  of  paper  width  and  places 
the  time  ordinate  down  the  paper.  As  many  pages  of  paper  as 
may  be  required  by  the  length  of  the  maneuver  are  then  used 
at  the  rate  of  66  time  points  per  page.  If  no  plots  are 
desired,  then  no  plot  data  cards  are  used. 

The  other  form  of  the  maneuver  data  is  the  page  of  data  which 
may  be  printed  at  every  time  point  or  at  every  mth  time  point. 

A.  sample  of  these  data  are  given  in  Table  IV.  These  data  have 
ueen  selected  as  being  most  pertinent  to  most  maneuvers.  The 
force  and  moment  summary  at  the  bottom  of  the  page  is  quite 
useful  in  determining  what  components  of  the  aircraft  are 


forces 

moments 

angles  of  attack 
aerodynamic 
coefficients 
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generating  the  forces  and  moments  acting  on  the  aircraft.  It 
should  be  remembered  by  the  user  that  the  rotor  forces  are 
strongly  affected  by  azimuth  position  and  will  vary  consider¬ 
ably  from  time  point  to  time  point  during  the  maneuver. 

3 .  Capabilities 

There  are  three  primary  types  of  inputs  which  may  be  simulated 
to  generate  maneuvers.  These  are:  control  inputs  by  the  pilot, 
control  inputs  by  simulated  automatic  systems,  and  external 
disturbances.  The  detailed  specification  of  these  inputs  are 
given  in  Volume  II,  Section  IIA  and  IiB.  Basically,  the  time 
histories  of  control  motion  or  external  disturbances  are  prede¬ 
termined.  The  operating  constants,  such  as  dead-band,  maximum 
and  minimum  rates,  thresholds,  etc.  of  automatic  systems  are 
also  predetermined  but  the  actual  control  output  of  these 
systems  depend  upon  the  values  of  position,  velocity,  and 
acceleration  during  the  maneuver.  Once  these  input  data  are 
defined,  the  resulting  aircraft  motion  is  computed  by  the 
maneuver  section  of  the  computer  program.  Any  combination 
of  the  above  inputs  may  be  operating  simultaneously  during 
the  maneuver.  If  certain  modes  of  motion  are  not  of  interest 
in  a  given  problem,  rotational  degrees  of  freedom  may  be 
effectively  locked  out  by  using  extremely  large  inertias 
about  the  appropriate  axes. 

There  are  nine  inputs  that  relate  to  pilot  control  actions. 
These  controls  are: 

1.  The  four  primary  controls 

(a)  Collective  pitch  lever 

(b)  Longitudinal  cyclic  stick 

(c)  Lateral  cyclic  stick 

(d)  Pedals 

2.  Power  controls 

(a)  Rotor  torque 

(b)  Auxiliary  thrust 

3.  Conversion  controls 

(a)  Change  mast  tilt  angle  (and  RPM) 

(b)  Change  tail  rotor  into  pusher  propeller 

(d)  Engage  rotor  brake. 
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There  are  eight  inputs  that  relate  to  functions  that  would 
normally  be  performed  by  automatic  systems.  These  are: 

1.  Rotor  speed 

(a)  Collective  pitch  governor 

(b)  Set  speed  for  governor 

2.  Rotor  attitude 

(a)  Flat  tracker  for  main  rotor 

(b)  Flat  tracker  for  tail  rotor 

(c)  Rotor  fold  mechanism 

3.  Autopilot  functions 

(a)  Yaw  stabilization 

(b)  Pitch  stabilization 

(c)  Sinusoidal  control  motion 

There  are  five  inputs  that  relate  to  external  disturbances. 

1.  Gusts 

(a)  Vertical  ramp-shaped  or  step-shaped  gust 

(b)  Vertical  sine-squared-shaped  gust 

(c)  Horizontal  ramp-shaped  or  step-shaped  gust 

(d)  Horizontal  sine-squared-shaped  gust 

2.  Weapon  recoil 
D.  ROTOR  AIRLOADS 

The  present  computer  program  was  derived  from  another  program 
(Reference  5)  which  contains  a  detailed  aeroelastic  represen¬ 
tation  of  the  rotor  blade.  This  portion  of  the  program  had 
to  be  deleted  because  the  mathematical  model  of  the  stop-fold- 
rotor  is  time-variant  and  not  quasi-static  like  the  model  of 
Reference  5.  This  time-variant  analysis  is  not  compatible 
with  the  dynamical  representation  of  the  rotor  used  in 
Reference  5.  The  airload  data  may  be  printed  out  at  up  to 
14  specified  time  points  during  a  maneuver.  These  data  are 
based  upon  rigid  blades  which  may  assume  different  coning 
angles  from  the  flapping  hinge  point.  No  provision  is  made 
for  the  effect  on  blade  loads  of  the  type  of  hub,  mass  distri¬ 
bution  of  the  blades,  or  aeroelastic  feedback.  Provisions 
for  possible  addition  of  these  effects  have  been  left  in  the 
program . 
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The  aerodynamic  parameters  of  usual  interest  are  printed  at 
each  of  20  blade  stations  for  twelve  equally-spaced  azimuth 
locations  in  the  trim  and  stability  sections  of the  program. 

The  maneuver  section  of  the  program  computes  parameters  at  b 
azimuth  stations,  where  b  is  the  number  of  blades.  Samples 
of  these  data  are  given  in  Table  V  for  a  case  with  no  radial 
flow  represented.  If  radial  flow  is  used,  then  two  lines  per 
blade  station  are  used  to  print  the  data  (see  Figure  35, 

Volume  II).  Complete  descriptions  of  these  data  are  given  in 
Volume  II,  Section  IIIJ. 

If  airloads  are  desired  only  for  the  trim  condition,  then  the 
maneuver  section  of  the  program  must  be  entered  and  airloads 
requested  at  time  t  =  0  seconds.  The  maneuver  should  be  set 
up  with  no  inputs  and  terminate  at  At  seconds.  If  airloads 
are  required  during  a  maneuver,  the  times  must  be  specified  in 
advance.  This  usually  requires  the  user  to  compute  the  maneu¬ 
ver  first,  examine  the  data  for  times  when  minimum  and  maximum 
values  of  pertinent  parameters  occur,  and  then  select  those 
times  for  airload  analysis  and  rerun  the  maneuver. 

E.  VECTOR  ANALYSIS 

The  primary  purpose  of  having  the  capability  to  oscillate  any 
one  of  the  pilot's  controls  in  a  sinusoidal  manner  is  to  enable 
the  response  computed  in  the  maneuver  section  to  be  reduced  to 
vector  form.  The  advantage  of  this  capability  becomes  apparent 
when  trying  to  isolate  cause  and  effect  relationships  of  new 
configurations  or  in  correlating  with  the  results  of  the 
linear,  uncoupled  stability  analysis. 

The  user  must  specify  at  what  time  point  the  vector  analysis 
is  to  begin  and  to  end.  Generally,  two  cycles  of  control 
excitation  is  needed  before  starting  the  analysis  and  the 
analysis  is  stopped  at  the  end  of  the  third  cycle.  These 
start-stop  times  may  vary  depending  upon  the  frequency  char¬ 
acteristics  of  the  aircraft  and  of  the  excitation.  The 
further  removed  the  excitation  frequency  is  from  a  natural 
frequency  of  the  aircraft,  the  sooner  the  vector  analysis  may 
start.  Any  parameter  that  may  be  plotted  may  also  be  vecto- 
rially  analyzed.  Basically,  all  that  is  done  is  to  fit  the 
computed  time  history  of  any  parameter,  P,  with  the  following 
equation.  The  coefficients  are  chosen  to  minimize  the  root- 
mean  -square  value  of  the  deviations  from  the  actual  curve. 

The  frequence,  w,  used  is  the  excitation  frequency  specified 
by  the  user. 

P  =  P  +  P.  cos  (tot  +  <t>  ) 
oi  p 

The  variance  is  also  computed  and  printed  out  for  each  para¬ 
meter  to  indicate  how  good  the  curve  fit  is.  The  variance 
should  be  between  0.95  and  1.00,  the  latter  of  which  repre¬ 
sents  a  perfect  fit.  A  sample  analysis  is  given  in  Table  VI. 
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Once  the  maneuver  parameters  have  been  determined  in  vector 
form,  each  may  be  compared  to  the  others  in  terms  of  phase 
angles  and  magnitudes.  This  enables  direct  comps rison  with 
the  same  data  computed  in  the  stability  section.  It  also 
indicates  the  degree  of  nonlinearity  and  coupling  that  exists 
for  the  specified  flight  condition  of  the  aircraft.  If  good 
agreement  exists  between  the  phase  angles  and  magnitudes 
computed  by  the  stability  and  the  maneuver  section,  then  the 
faster,  linear,  uncoupled  analysis  of  the  stability  section 
may  be  assumed  to  be  representative  of  the  characteristics  of 
the  aircraft. 

Another  useful  relationship  may  be  determined  with  the  para¬ 
meters  in  vector  form.  Any  parameter  may  be  expressed  as  a 
linear  combination  of  two  other  parameters  provided  that 
neither  of  the  three  parameters  differ  by  0  or  180°  in  phase 
angle.  This  feature  is  often  useful  in  isolating  cause  and 
effect  relationships  of  different  phenomena  exhibited  by  new 
configurations.  An  example  of  this  analysis  is  given  in 
Table  VI. 
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SECTION  III 

DESCRIPTION  OF  THE  MATHEMATICAL  MODEL 

This  section  describes  the  important  equations  that  are  used 
in  the  mathematical  model.  It  is  not  intended  to  be  an  ex¬ 
haustive  treatment;  the  FORTRAN  listing  in  Volume  IV  is  the 
controlling  document. 

A.  FUSELAGE  AERODYNAMICS 

All  aerodynamic  data  of  the  fuselage  are  represented  in  the 
stability-axis  system.  If  wind-tunnel  data  on  the  rotorcraft 
to  be  simulated  are  available,  the  force  and  moment  data 
should  be  resolved  to  stability  axis.  The  aerodynamic  center 
of  the  fuselage  must  be  specified.  A  sample  of  such  data  from 
a  wind-tunnel  test  is  plotted  in  Figure  3.  The  lines  indicate 
the  fairing  for  determining  the  coefficients  used  in  the  input- 
data  cards  of  the  computer  program.  Since  no  provisions  are 
made  for  rolling -moment  data,  the  vertical  distance  between  the 
center  of  gravity  and  the  aerodynamic  center  of  the  fuselage  is 
determined  so  that  the  fuselage  force  produces  the  required 
rolling  moment.  Other  fuselage  data  are  input  directly  from 
the  fitted  curves  as  indicated  in  Figure  3.  If  curves  are 
nonlinear,  the  slope  near  the  point  of  interest  should  be  used. 
For  example,  if  pedal  control  margins  are  being  investigated 
the  curve  should  be  approximated  at  high  sideslip  angles  in¬ 
stead  of  zero.  If  level  flight  stability  data  were  of  inter¬ 
est,  then  the  curves  should  be  fitted  near  zero  sideslip.  If 
wind-tunnel  data  are  not  available,  the  aerodynamic  center  and 
other  aerodynamic  data  must  be  estimated,  giving  proper  con¬ 
sideration  to  the  representation  of  forces  and  moments. 

The  mathematical  model  of  the  t ilting-proprotor  aircraft  in¬ 
cludes  a  representation  of  the  movement  of  the  center  of  gravity 
as  the  proprotors  convert  from  helicopter  to  airplane  mode. 

The  change  in  drag  area  and  the  pitching  moment  induced  by 
pylon  drag  are  also  represented.  The  input  value  of  center-of- 
gravity  location  assumes  that  the  pylons  are  in  the  helicopter 
configuration.  The  value  of  pylon  flat-plate  drag  area  re¬ 
presents  the  difference  in  drag  with  the  masts  vertical  and 
horizontal . 

The  equations  used  in  the  representation  of  fuselage  charac¬ 
teristics  are  summarized  in  Table  VII. 

To  simplify  the  nomenclature,  the  variables  that  are  taken 
from  the  input  data  are  represented  by  their  FORTRAN  names. 

The  user’s  guide  (Volume  II)  defines  these  symbols  in  more 
detail . 
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B.  AIRFOIL  AERODYNAMICS 

The  standardized  input  data  on  three  cards  are  used  to  compute 
the  coefficients  of  lift,  Cl,  and  drag,  Cq,  as  a  function  of 
angle  of  attack  and  Mach  number.  In  the  following  discussion 
YXXI  refers  to  the  Ith  aerodynamic  input  for  the  appropriate 
aerodynamic  surface  YXX,  as  defined  in  Table  VIII.  Except  for 
sections  specifically  labeled  otherwise,  formulas  and  proce¬ 
dures  apply  to  all  aerodynamic  surfaces;  i.e.,  rotors,  wings, 
elevators,  or  fins. 

I.  Lift  Characteristics 

If  the  lift  slope  for  zero  Mach  number,  input  YXXI7,  is  zero, 
the  program  will  compute  Cl  and  Cd  from  data  tables  for  a 
64A210  airfoil.  If  input  YXXI7  has  a  nonzero  value,  Cl  and 
Cd  will  be  computed  from  the  analytical  and  empirical  equa¬ 
tions  described  in  the  following  discussion.  The  airfoil  sec¬ 
tion  characteristics  included  in  the  aerodynamic  input  listed 
in  Table  VIII  are  for  a  two-dimensional  symmetrical  section  at 
zero  Mach  number.  These  input  values  are  corrected  by  the 
program  to  obtain  the  three-dimensional  characteristics  for 
an  unswept,  symmetrical  surface.  The  form  of  the  computed 
three-dimensional  characteristics  are  illustrated  in  Figure  4. 

a.  Angle  of  Attack  Effects 

The  current  equations  represent  only  symmetric  airfoil  sections. 
The  same  equations  may  be  used  for  positive  and  negative  angles 
of  attack;  since  for  symmetric  sections  negative  angles  of 
attack  produce  the  same  magnitude  of  lift  coefficients  as 
positive  angles,  only  of  opposite  sign.  Thus,  only  angles  of 
attack  between  0  and  180  degrees  need  be  considered  when 
evaluating  the  operations  for  Cl.  This  angle  of  attack  range 
is  divided  into  two  major  flow  regions,  normal  (0°  <  a  <  90°) 
and  reverse  (90°  <  a  <  180°)  flow,  which  utilize  the  same 
equations  but  require  separate  input  data. 

The  normal  and  reverse  flow  regions  are  further  divided  into 
three  operating  regions: 

(1)  The  region  where  lift  is  linearly  proportional 
to  angle  of  attack,  i.e.,  0°  <  a  <  astall* 

(2)  The  stalled  region  where  flow  separation  has 
occurred  over  part  of  the  upper  surface,  i.e., 
astall  <  a  ^  as. 

(3)  The  region  where  the  surface  can  be  considered 
to  be  a  flat  plate,  i.e.,  cig  <  a  <  90°. 
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The  boundaries  of  these  regions  defined  as  a^tiall  an<*  aB  are 
determined  internally  in  the  program.  The  first  of  these, 
astall»  is  defined  as  the  three-dimensional  angle  of  attack  at 
which  Cl  =  CLmax  for  the  Mach  number  under  consideration.  The 
other  limit,  ag,  is  defined  as  the  stall  angle  of  attack  at 
M  =  0  plus  5  degrees.  If  ag  is  greater  than  40  degrees,  it  is 
reset  to  40  degrees,  and  CLmax  at  M  =  0  is  adjusted  so  that 
astall  =  35  degrees. 

b.  Mach  Number  Effects 

To  provide  the  best  approximation  of  experimental  data,  dif¬ 
ferent  sets  of  equations  are  used  to  compute  Cl  in  the  sub¬ 
sonic,  transonic,  and  supersonic  regions.  The  boundaries 
of  three  Mach  regions  are  defined  as: 

(1)  Subsonic:  0  <M  <  YXX1 

(2)  Transonic:  bounded  by  YXX1  on  the  lower 
end,  and  the  maximum  of  either  YXX2  or  the 
value  of  Mach  number  that  satisfies  the 
following  condition 

(YXX18)  ^M^  -  1  =  1  on  the  upper  end. 

(3)  Supersonic:  bounded  only  on  the  lower  end 
by  the  value  coincident  with  the  upper 
boundary  of  the  transonic  region. 

The  value  of  Cl  below  the  stall  angle  is  computed  as  a  linear 
function  of  the  two-dimensional  lift  curve  slope,  a2p,  and 
the  effective  angle  of  attack,  a, 

CL  =  a2Da 

where  a  =  aQ  -  and  oq  is  the  geometric  angle  of  attack  and 
a£  is  the  induced  angle  of  attack. 

The  program  computes  the  lift  curve  slope  by  one  of  the 
following  equations,  depending  upon  the  Mach  number.  The 
two-dimensional  lift  curve  slope,  YXX17,  is  corrected  for 
compressibility  effects  at  subsonic  speeds  by  the  Prandtl- 
Glauert  correction: 

_  YXX17 

2Dc '  7rrw 

For  the  supersonic  region,  the  lift  curve  slope  is  computed 
from  the  small  perturbation  theory  equation: 

_ _ 1_ 

2D 
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For  the  transonic  region,  the  value  of  lift  curve  slope  is 
determined  from  a  second-order  polynomial  curve  in  Mach  number 
that  satisfies  the  following  conditions: 

(1)  The  value  of  the  lift  curve  slope  at  the 
boundaries  of  the  region  must  equal  the  values 
computed  at  these  points  by  the  equations  for 
the  adjacent  region. 

(2)  The  slope  of  the  lift  curve  slope  versus  Mach 
number  curve  at  the  common  boundary  of  the 
transonic  and  supersonic  region  must  be  continuous 
with  Mach  number. 


These  constraints  on  the  lift  curve  slope  provide  continuous 
values  of  lift  curve  slope  with  Mach  number,  and  result  in  a 
discontinuity  in  the  slope  of  the  lift  curve  slope  versus  Mach 
number  curve  only  at  the  input  Mach  number  of  YXX1. 


The  lift  coefficient  for  the  angle  of  attack  region  where  the 
surface  can  be  considered  to  be  a  flat  plate,  ag  <  a  <  90°, 
is  computed  from  the  following  empirically  derived  equation: 


[( 


2K^sina 


-  K 


2)k3 


+  0.81  cosa 


The  values  of  the  constants  and  K2  are  selected  internally, 
with  different  values  being  set  for  the  rotor  and  for  fixed 
surfaces. 


0.938  for  the  rotor 

YXX3  or  YXX7  for  fixed  surfaces  for 
normal  and  reverse  flow  respectively 

0.581  for  the  rotor 

0.81  for  fixed  surfaces 


The  constant  K3  is  selected  independently  of  the  type  of  sur¬ 
face,  but  is  computed  from  one  of  the  following  equations 
depending  upon  the  Mach  number  and  angle  of  attack. 

1  +  O^SM^  if  M  £  1  and  a  >  ag  or  for 

computing  the  value  of  ag 

0 . 84  +  M°-j-Q-Q-Q  if  M  >  1  and  a  >  aB 


The  lift  coefficient  in  the  angle  of  attack  range  between 
astall  and  aB  is  determined  by  linear  interpolation  between 
the  values  of  Cl  at  astall  and  Cl  at  ag. 
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c .  Induced  Angle  of  Attack 

In  all  of  the  Mach  number  and  Angle  of  attack  regions,  the 
induced  angle  of  attack  is  considered  to  be  proportional  to 
the  coefficient  of  lift.  The  expression  for  induced  angle  of 
attack  was  derived  from  PrandtI's  momentum  theory  as  modified 
to  better  represent  experimental  data.  This  expression  is  in 
the  form  of  one  of  the  following  equations: 


„  _  /  YXXI7  \ 

ai  ~  - 2 -  a 

1  \YXXI7  +  4tiYXXI8/ 


if  a  <  a 


stall 


a  .  = 


1  jti(YXXI8) 


if  a  >  a 


stall 


For  elevators  and  fins  an  additional  correction  is  made  to  cL 
to  account  for  the  effects  of  tail  boom  bending.  This  correc¬ 
tion  is  accomplished  by  the  following  equation: 


i  •  nxxu^jjvi, 


where 

At  =  area  of  the  surface 

qt  =  dynamic  pressure  at  the  surface. 

2 .  Drag  Characteristips 

The  equation  utilized  for  the  computation  of  the  coefficient 
of  drag,  Cq,  is  dependent  on  the  Mach  number.  The  three  Mach 
ranges  for  Cp  are:  subsonic,  Os  M  <  I;  transonic,  1  <  M  <  Ms; 
and  supersonic,  M  >  Ms,  where  Ms  is  the  lower  boundary  of 
the  supersonic  range. 

a .  Unstalled  Flow 

At  subsonic  Mach  numbers,  the  airfoil  is  considered  to  be 
operating  either  unstalled,  0  <  a  <  asx,  or  as  a  stalled  flat 
plate,  asx  <  a  <  90°,  where  asx  is  defined  as  the  minimum  of 
either  astaIL  or  the  value  of  a  that  satisfies  the  condition 
that  the  unstalLed  drag  equation  be  equal  to  the  input  YXXL6. 
This  second  limiting  condition  is  applied  so  that  the  subsonic 
value  of  Cj)  does  not  exceed  the  transonic  value  of  Cp  and 
result  in  a  discontinuity. 
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When  0  <  a  <  asx,  the  drag  coefficient  is  computed  from  a 
second  order  polynominal  in  a,  the  coefficients  of  which  are 
input  data.  The  equation  used  is: 

C„  =  -Y-— J-.?  +  YXX13a  +  YXX14a2  +  CD3 

D 

The  term  CD3  is  an  empirically  derived  term  added  to  Cp  to 
account  for  the  effects  of  Mach  number  on  Cp.  This  term 
results  in  the  characteristic  drag  divergence  with  increasing 
angle  of  attack  and  Mach  number.  This  term  is  defined  as: 


CD3 


max 


0.0 


max 


M 

0.35 


YXX1  +  0.0332a 


b.  Stalled  Flow 


When  asx  <  a  <90,  it  is  necessary  to  use  separate  equations 
for  the  rotor  and  for  fixed  airfoils.  For  fixed  airfoils  such 
as  the  elevators  and  fins,  Cp  may  be  represented  by: 


CDsx  -  2.0 

7a  -  90) 2 
sx  J 


+  2.0 


This  empirical  equation  yields  a  parabolic  curve,  which  is 
concave  downward.  It  has  a  value  of  CDSX  at  as^  and  a  maximum 
value  of  2.0  at  90  degrees.  The  value  of  CDSX  is  the  minimum 
of  either  YXX16  or  the  value  of  Cp  computed  from  the  unstalled 
equation  for  a§x.  The  maximum  of  Cp  at  90  degrees,  2.0,  was 
selected  as  being  representative  of  a  flat  plate  perpendicular 
to  the  flow.  Because  of  the  limiting  value  for  at  90 
degrees,  it  is  necessary  to  limit  the  value  of  YXX16  to  values 
less  than  or  equal  to  2.0  in  order  to  maintain  the  proper 
shape  of  the  CD  versus  a  curve. 


For  the  rotor  sections  operating  in  the  a  range  where 
agx  <  a  <  90° ,  the  equation  for  CD  is: 


where 


CD  =  CD4sin 


i.n2a  +  £ 


CD  -  CD4sin2 
sx 


cosa 


sx 


oleosa 


CD4 


2.1 


1  +  0.25M4  for  M  £  1 

°-81* +  m0:T8  for  M  > 1 
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For  the  range  of  Mach  numbers  where  M  >  1  but  less  than  the 
lower  boundary  of  the  supersonic  region,  Cp  has  a  constant 
value,  equal  to  YXX16,  regardless  of  the  value  of  a. 

When  the  value  of  Mach  number  is  greater  than  the  lower 
boundary  of  the  supersonic  region,  the  value  of  Cj),  for  all 
values  of  a,  is  the  minimum  of  either  YXX16  or  that  determined 
by  the  following  equation: 


C 


_  4(a  +  YXX15) 

D  "  /M2~-T 


+  YXX12 


C .  ROTOR  REPRESENTATI  ON 

The  aerodynamic  characteristics  of  two  rotors  may  be  repre¬ 
sented  in  the  computer  program.  Depending  on  the  input-data 
values,  the  two  rotors  can  be  used  as 


A  vertical-thrusting  main  rotor  and  a  right- 
thrusting  tail  rotor 

Side-by-side,  contra-rotating  main  rotors 
Tandem,  contra-rotating  main  rotors 


In  the  latter  two  configurations,  the  input  data  for  the  Tail- 
Rotor  Group  applies  to  the  left  or  aft  rotor.  Much  of  the 
input  data  and  all  the  analysis  are  the  same  for  both  rotors, 
so  the  same  description  applies  to  both.  All  the  aerodynamic 
equations  are  written  in  the  mast-axis  system.  Infinitely 
stiff  blades  are  assumed,  with  a  variety  of  available  hub 
geometries : 


-  Rigid,  with  a  virtual  hinge  point 

-  Spring-restraint  (nonisotropic) 

-  Gimbaled 

-  Flapping-hinge-offset 

-  Focused-pylon 
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The  basic  method  used  to  compute  the  rotor  forces,  moments, 
and  power  required  uses  conventional  blade-element  theory,  w  »8) 
A  maximum  of  20  evenly-spaced  blade  stations  may  be  repre¬ 
sented.  Twelve  azimuth  stations  are  used  in  the  trim  and 
stability  sections  of  the  program.  The  maneuver  section  uses 
as  many  azimuth  stations  as  blades.  The  number  of  blades  (of 
rectangular  planform),  as  well  as  the  blade  radius  and  chord, 
are  specified. 


1.  Special  Provisions 

Each  blade  segment  is  assumed  to  be  untwisted,  but  each  may 
have  a  different  incidence  with  respect  to  the  hub.  Linear 
twist  is  represented  by  20  finite  values  of  blade  pitch  which 
are  input  on  the  fifth,  sixth  and  seventh  data  cards  in  the 
rotor  group.  The  values  to  be  input  are  given  by: 


0 


0 


root 


♦  e(x) 


where  the  inputs  are  made  from  the  tip  of  the  blade  (x  =  1)  to 
the  root  (x  =  0).  Since  the  pitch  of  most  rotors  is  less  at 
the  tip  than  at  the  root,  the  input  values  are  normally 
negative.  For  example,  the  UH-1B  rotor  has  a  10-degree 
"washout"  of  pitch  from  root  to  tip.  Thus,  for  a  20-segment 
rotor , 


XMR29  = 

-10.0 

(blade  tip,  x 

XMR30  = 

-  9.5 

(x  =  0.95) 

XMR31  = 

• 

• 

-  9.0 

(x  =  0.90) 

XMR47  = 

-  1.0 

(x  =  0.10) 

XMR48  = 

-  0.5 

(x  =  0.05) 

Precone  is  accurately  represented  by  the  program.  The  input 
value  is  the  amount  of  precone  built  into  the  hub.  The  actual 
coning  angle  is  computed  using  a  quasi-static  balance  among 
centrifugal-force  moment,  airload  moment,  and  coning-spring 
moment,  represented  as  follows: 


n 


2 

i=l 


rbL 


H 


n 


S^sinp 


«£  [(rb)  C0spo +  erb]ral(irb)  ‘ 


KcX0 


=  0 


i=l 


The  Newton-Raphson  iteration  technique  is  used  to  find  a  new 
value  of  Xq  for  each  pass  through  the  rotor  analysis.  After 
convergence,  Xq  represents  the  coning  deflection  of  the  blade 
from  that  value  of  coning  angle  which  is  specified  by  the 
input  data  thus,  p0  =  PpCinput  +  x0* 
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Provisions  are  made  for  blade-pitch  change  with  coning  as  well 
as  with  flapping  relative  to  the  mast.  The  use  of  hub-drag 
coefficient  improves  the  correlation  of  computed  and  flight- 
test  data  for  stick  positions  and  gradients  in  forward  flight. 

A  first-order  approximation  of  a  nonuniform  induced-velocity 
distribution  is  possible  by  use  of  a  coefficient  for  tip- 
vortex  effect.  This  effect  does  not  influence  total  rotor 
forces  and  power,  but  it  does  affect  the  spanwise  and  azimuthal 
distribution  of  airloads.  It  is  represented  by: 


V.  =  Jv  2  +  V  2  (velocity  of  the  hub  in  the  X-Y  plane) 

■  ■  a  y 

Depending  on  the  values  of  and  S2,  a  value  of  XK  is  deter¬ 
mined  for  later  use  by  whichever  equation  below  is  true  first. 

XK  =  0.5  if  52  <  1 


XK  =  11.25  Vh/S2R  if  Vh/S2R  <  0.1067 

XK  =  1.36  -  1.5  Vh/52R  if  Vh/52R  <  0.5733 

XK  =  0.5  if  Vh/S2R  *  0.5733 

A  correction  factor  to  the  momentum  theory  value  of  induced 
velocity  is  then  computed  as  follows: 


xy  =  V.  XK(XMR27)'\/-  'C  +  \/ £  +  Vj2 

mom  V  z  V 

Rotor  induced  velocity,  Vi_0_,  is  determined  from  the  induced 
velocity  subroutine  (9)  which  is  based  on  momentum  theory  (10) 
with  empirical  corrections  so  that  a  value  of  induced  velocity 
is  assured  for  all  rotor  operating  states. 

The  input  variable  XMR27  should  be  zero  if  tip  vortex  effect 
is  not  desired  and  should  be  unity  if  it  is  desired.  The 
second  term  of  the  equation  below  which  involves  xy  is  imposed 
upon  the  induced  velocity  distribution  only  in  the  azimuth 
intervals  between  45  to  105  degrees  and  255  to  315  degrees  and 
only  outboard  of  the  70%  blade  station. 

V.  =  V.  |"l  -  S*(l  ♦  XKcos\J/)l  ♦  xy  sin(6(\J/  -  45°)] 
mom L  J 

where  x  is  the  fraction  of  blade  radius,  rb/R. 

2  .  Blade -Element  Aerodynamics 

The  angle  of  attack,  a,  and  the  sweep,  y,  at  a  specified  blade 
station,  r,  and  azimuth  station,  \J/ ,  is  computed  using  the 
following  equations: 


23 


a  =  0  +  4) 


where 

0  =  Aq  -  tan“^cos\|/tanA^  +  sin\J/tanB^J  -  ^p-p^JtanS^ 

and  4>  is  determined  by  resolving  wind  vectors  as  follows: 

The  three  components  of  velocity  at  a  blade  element,  illus 
trated  in  Figure  5,  are  computed  by: 


Up  =  XS2r,  cosp  -  r.(?  -  (V  cos\J /  -  V  sin\J')sinp  +  UD 

D  D  xh  yh  gust 

UT  =  S2r,  cosp  ♦  V  sinvl/  ♦  V  cos\J/  -  UT 
10  xh  yh  Agust 

UD  =  -AS2r,  sinp  -  (V  cos\j/  -  V  sin\J/)cosP  +  UD 
R  b  xh  yh  Rgust 

where  X  =  (V  -  V.)/S2R 
zh  1 

Gust  velocities  are  determined  by  the  following  equations: 

Letting  Gjj,  Gy,  and  Gy  b  the  longitudinal,  vertical,  and 
lateral  velocities  of  th»  gust,  respectively,  and  defining 
GI  =  ghcos^  “  Gys^n^»  the  Sust  components  are  given  by: 


UP 

=  GjSinp  -  Gycosp 

gust 

UT 

=  Gjjsin\J/  -  GyCosty 

gust 

UR 

=  Gjcosp  +  Gysinp 

gust 

The  variables  Vx^,  Vy^,  and  Vz^  are  the  velocities  at  the  rotor 

hub  caused  by  airframe  linear  and  angular  velocities,  mast  ori¬ 
entation,  and  aircraft  attitude  in  space.  Let  the  components  in 
body  reference  of  free-stream  velocity  be  Vx,  Vy,  and  Vz;  the 
angular  rates  of  the  body  be  p,  q,  and  r;  and  the  linear  dis¬ 
tances  from  the  center  of  gravity  to  the  rotor  hub  be  lx^,  iy^* 
and  lzh,  all  of  which  are  referenced  to  the  X,  Y  and  Z  1 

axes,  respectively.  The  linear  velocity  at  the  rotor  hub 
caused  by  angular  velocities  is  given  by 
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V  =  ql  -  r 1 
xr  zh 


Vy  =  rlx.  -  PLZ, 
7r  h  1 


=  pi..  -  qi. 


To  these  values  are  added  the  free  stream  velocity  components. 
The  total  velocity  in  body  reference  is  then: 


=  V 

xr 

♦  V 

xbody 

=  V 

♦  V 

^body 

=  V 

+  V 

zr 

zbody 

where  VyK  ,  is  corrected  for  the  input  value  of  sidewash  co¬ 
efficient  for  the  tail  rotor  of  the  conventional  helicopter 
configuration. 

The  final  transformation  requires  the  orientation  of  the  mast 
relative  to  the  body  to  be  known.  Mast  yaw  angle  is  currently 
set  to  zero.  Mast  pitch  angle  is  determined  by  input  data  as 
pretilt,  for  a  conventional  helicopter,  or  conversion  angle,  for 
the  tilt-proprotor  aircraft,  and  is  zero  when  perpendicular  to 
the  X  -  Y  plane  of  the  aircraft.  Mast  roll  angle  is  zero  for 
the  tilt-proprotor  and  tandem  configurations  and  for  the  main 
rotor  of  the  conventional  helicopter  configuration.  The  roll 
angle  is  90  degrees  for  the  tail  rotor  of  the  conventipnal 
helicopter.  Defining  and  ©m  as  the  mast  roll  and  pitch 
angle,  the  free  stream  velocity  components  at  the  hub  in  the 
shaft  reference  system  are  given  by  the  following  equations. 

V  =  V  cosG  -  V  sine 

x  m  z  m 

V  =  V  sin<t>  sin0m  ♦  V  „cos<t>  +  V  cos0  sin<t> 
y,  xmmy  m  zmm 

V„  =  V  cos<{>  sin0  -  V  sin4>  +  V  cos©  cos4> 
zh  x  m  m  y  m  z  m  m 
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These  velocities  are  then  used  in  the  preceding  equations  for 
Up,  UT,  and  Ur  at  the  blade  element.  The  inplane  and  total 
velocities  are  then  calculated  by: 

UI  -  ^  UT2  +  UR2 
u  =  /uj2  +  Up2 

Radial  flow  is  treated  by  obtaining  two  lift  and  two  drag 
coefficients.  One  pair  of  coefficients  is  resolved  in  the 
conventional  manner  for  rotor  analysis.  The  other  pair  is 
associated  with  radial  flow  and  act  as  if  the  freestream 
velocity  were  parallel  to  the  blade.  The  conventional  angle 
of  attack  ic,  determined  using  an  inflow  angle,  <J>,  from  the 
following  equation.  Defining  sweep  angle,  y,  as: 


tany  = 


the  inflow  angle  is 
and  so , 


Up  o 

tan<J>  =  rr—sin 
UT 


r 


a  =  e  -  0 


M 


U  .  2 
=  ^sin 


r 


where  <J>  is  placed  in  the  proper  quadrant ,  depending  upon  the 
sign  of  Up  and  Ur.  For  no  radial  flow,  Ur  is  zero  and  y  is 
90  degrees  so  the  equation  for  4>  is  the  same  as  conventional 
rotor  theory.  The  lift  coefficient,  C^,  and  drag  coefficient, 
Cp,  are  obtained  from  either  the  airfoil  data  tables  or  from 
tne  equation  subroutine  described  in  the  preceding  section  of 
this  report. 


The  lift  and  drag  coefficients  associated  with  radial  flow  are 
obtained  from  the  same  data  tables  or  equations.  The  angle  of 
attack  for  radial  flow  is  defined  as: 


tana^ 


Up  2 

rs-cos  y 
UR 


and  the  Mach  number  is 


After  the  lift  coefficient  is  obtained  it  is  modified  by  the 
input  value  in  XMR26,  which  should  be  between  0  and  10  as 
follows 
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„  _  „  1”  XMR26  2*  c  1 

clr  ~  l  [37.3(c  V  4r>J 

where  c  is  the  rotor  chord  length. 

The  radial  lift  and  drag  coefficients  are  represented  by  Clr 
and  Cj)R  respectively.  All  lift  and  drag  vectors  are  assumea 
to  be  aligned  perpendicular  and  parallel  to  the  local  relative 
wind. 

The  lift  and  drag  coefficient  vectors  are  then  resolved 
parallel  and  perpendicular  to  the  blade  to  obtain  net  coeffi¬ 
cient  vectors  in  the  lift,  drag,  and  radial  directions.  The 
equations  used  for  the  nth  blade  segment  are  as  follows : 


= 


= 


A„  = 


=  U 


.  2 
sin 


r  (cLUT  ♦  CDUpsin2r)  ♦  cos2rcose/cL  UR  +  CD  Upcos2yj 

\  R  R  - 


U(CDUT  "  CLUpsin2r)sin2Y 

Ucos2ycosG [C^  UD  -  CT  U0cos2y) 
\  °R  R  LR  P  / 


Total  rotor  forces  and  power  are  obtained  by  summing  these 
airloads  along  each  blade  as  follows: 

N.  N 

\|/  r 


Thrust  =  £  £(ALcosp  -  ^sinf) 


r  *  n=l 
N*  Nr 


H-Force  =  jfr  ^  ^  jApsinil'  -  ^ALsinp  ♦  ARcospJ  cos\!' j 

r  *  i=i  n=lL 


N*  Nr 


Y-Force  =  ^  ^ j-  ApCosi)  -  (ALsinp  +  ARcospJ  sinil'  J 

r  ^  i=l  n=l 


N*  Nr 


Torque  =  Inr  ][  2  ADr 


r  *  &  n=l 
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Nr  is  the  number  of  blade  stations  and  is  the  number  of 
azimuth  stations.  equals  twelve  in  the  trim  and  stability 

sections  of  the  program  and  is  equal  to  the  number  of  blades 
in  the  maneuver  section. 

3  .  Blade  Motion  Equation 

The  basic  differential  equation  of  blade  motion  considers  the 
effects  of  shaft  angular  rates,  blade  weight  and  inertia, 
nonisotropic  hub  spring  restraint,  offset  flapping  hinges,  and 
blade  airloads.  The  flapping  acceleration  vector  is  computed 
for  each  blade  from  the  following  equation. 


pj  =  (cosp  +  q)  ^pmcos\J/  -  qmsin\l/)cosp 


-Uk' 


lF/A 


cos  ♦  K.  sin 


^LAT 


in2\J/ 


(gyroscopics) 


(spring  restraint) 


-  S2z(cos(?  +  q)sin(3 


(centrifugal  force) 


+  -pcR_  rA  r 

Vr  ^  Li  bi 

Nr  r  / 

-  ^  wirb  .  s^nP  (  sin\|/cos9r 


(airload  moment) 


ground 


sinl 


lground 


+  cos\J/sin0 


+  cos(3cos0 


ground 


ground 


cos<t> 


ground 


(weight  moment) 


e  r 

where  rj  =  •=—  I  r.  mdr 

*b  J0  b 

The  blade  accelerations  are  added  vectorially  to  obtain  a 
resultant  acceleration  magnitude  and  orientation  of  the  rotor 
disk.  The  projections  of  this  resultant  vector  on  each  blade 
is  then  used  to  determine  the  flapping  acceleration  of  the 
blade.  This  acceleration  for  the  ith  blade  is  given  by: 


p\  =  ^  p\cos  ( j  -  i )  A\J/ 

J=1 

where  A\J/  is  the  azimuth  angle  between  blades. 
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4 .  Geometry 

The  location  of  the  rotor  hub  must  )  e  known  for  purposes  of 
resolving  forces  about  the  center  of  gravity.  This  requires 
the  mast  pivot  point  and  mast  length  to  be  specified  for  con¬ 
figurations  which  use  tilting  rotors.  For  conventional 
helicopter  configurations,  the  v.iaft  pivot  point  location 
should  be  the  rotor  hub  location  and  the  mast  length  should 
be  input  as  zero.  The  mast  tilt  angle  must  be  correctly 
specified  so  that  flapping  is  accurately  represented. 

The  blade  weight  distribution  needs  to  be  accurately  repre¬ 
sented  so  that  the  correct  center  of  gravity  shift  is  computed 
during  blade  folding.  Also,  the  total  blade  inertia  about  the 
flapping  axis  must  be  accurately  represented  if  meaningful 
rotor  stability  derivatives  with  respect  to  angular  rates  are 
to  be  obtained.  A  prior  version  of  the  computer  program  (1»4) 
included  a  dynamic  analysis  of  the  rotor  blade  and  included 
blade  deflections  and  deflection  velocities  in  the  aerodynamic 
representation.  This  capability  of  the  program  had  to  be 
deleted  in  order  to  accommodate  the  additions  required  to 
represent  the  stop-fold  rotor  concepts.  The  blades  are 
currently  represented  as  rigid  airfoils  that  are  free  to  flap 
but  not  deflect. 

D.  FIXED  AERODYNAMIC  SURFACES 

All  of  the  fixed  aerodynamic  surfaces  represented  in  the 
computer  program  utilize  lifting  line  theory  with  corrections 
made  for  aspect  ratio.  The  center  of  pressure  location  of 
each  surface  must  be  specified  by  the  input  data.  The  airfoil 
aerodynamics  subroutine  discussed  previously  is  used  to  deter¬ 
mine  lift  and  drag  coefficients  at  the  angle  of  attack  and 
Mach  number  of  the  airfoil.  The  wing  and  elevator  are  assumed 
to  be  mounted  horizontally  and  the  fin  is  assumed  to  be 
mounted  vertically. 

1 .  Wing 

The  wing  area,  aspect  ratio,  two-dimensional  lift  curve  slope, 
maximum  lift  coefficient,  angle  of  incidence  of  the  zero  lift 
line  relative  to  the  fuselage  waterline,  and  center  of  pres¬ 
sure  location  of  the  right  panel  are  the  more  important  input 
parameters.  Wing  side  force  and  pitching  moment  are  not 
represented  by  the  computer  program.  If  a  cambered  airfoil 
section  is  used,  the  resultant  pitching  moment  should  be 
included  in  the  data  of  the  fuselage. 

The  wing  area  used  should  be  the  total  area  of  the  wing  as 
shown  in  a  top  view,  including  fuselage  carry- through.  In 
the  event  that  relatively  "stubby”  wings  are  used  and  no  wind 
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tunnel  data  are  available  to  provide  corrective  data,  the 
following  correction  to  the  slope  of  the  wing  lift  curve 
should  provide  an  accurate  representation  of  first-order 
wing-body  interference  effects :( 11) 


where  bf  is  the  width  of  the  fuselage  and  is  the  span  of 
the  wing.  The  above  value  of  lift  curve  slope  should  be  used 
for  the  two-dimensional  lift  curve  slope  which  is  specified 
by  the  input  data.  If  no  wing  is  used,  a  zero  should  be  input 
as  the  value  of  wing  area. 

The  wing  contributions  to  roll  and  yaw  moment  are  computed 
using  the  values  of  stability  derivatives  which  are  input  to 
the  program.  If  these  stability  derivatives  are  not  known 
from  flight  or  wind  tunnel  test,  they  may  be  easily  estimated 
using  standard  techniques.  O-*, 13)  These  derivatives  may  be  set 
to  zero,  if  no  wing  is  used,  or  if  only  longitudinal  stability 
is  of  interest  for  a  winged  aircraft.  The  roll  and  yaw  moment 
of  the  wing  is  computed  by: 


t  = 


b 

w 

w 


w 


N. 


w 


=  ftA,  |^XWG12  +  XWG13  CL  +  (xWG14  CjT  +  XWG15  p)t*j 

=  £SwV2bw  J|xWG16  ♦  XWG17  CL2jpf  +  j  |^XWG18  CL2 
+  XWG19-g^jr  +  |VwG20  CL  ♦  XWG21^YWG12 


+  2  YWG13 


aw5B~7-j)-01745]p  t*' 

w  f  -j 


The  effect  of  the  rotor  downwash  on  the  wing  is  represented 
by  adding  the  induced  velocity  of  the  rotor  to  the  Z-  and  X- 
component  of  velocity  which  is  used  to  determine  the  wing 
angle  of  attack.  The  calculated  value  of  induced  velocity  is 
multiplied  by  the  number  input  on  the  data  cards  and  is 
assumed  to  act  parallel  to  the  mast. 


Wing  downwash  and  dynamic  pressure  reduction  in  the  wake  is 
represented  by  the  following  equations.  The  deflection  of 
the  wing  wake,  in  the  direction  opposite  to  lift  is  given  by: 
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e  =  XWG09  CT 

W 

W 


The  dynamic  pressure  reduction  along  the  wake  centerline  is: 


-  XMG10/°P - 2/nD 


’max 


cos 


VZK 


) 


where  D  is  the  distance  from  the  wing  trailing  edge  to  the 
elevator  leading  edge  in  wing  chords,  h  is  the  half-width  of 
the  wing  wake  at  a  point  £  wing  chords  down  the  wake  center- 
line  .  The  wingspan  and  chord  are  not  input  directly  but  are 
calculated  from  the  input  values  of  area  and  aspect  ratio  as 
follows: 

cw  =V^, 

b  =  MTS', 

W  v  w  w 

When  no  wing  is  used,  cw  and  b*,  are  set  to  urity  for  purposes 
of  defining  t*  which  is  used  to  nondimensionalize  time  in  the 
stability  section  of  the  program. 

2 .  Horizontal  Stabilizer 

The  geometric  and  aerodynamic  specifications  of  the  hori¬ 
zontal  stabilizer  are  similar  to  those  of  the  wing  although 
elevator  wake  calculations  are  not  made.  Four  coefficients 
are  currently  used  to  represent  the  effects  of  the  wing  and 
rotor  wake  on  the  elevator.  XEL11  is  used  to  multiply  the 
wing  lift  coefficient  and  the  product  represents  the  down- 
wash  angle  at  the  elevator  caused  by  wing  lift.  The  rotor 
induced  velocity  is  added  to  the  Z-component  of  elevator 
velocity,  after  multiplication  by  a  constant,  in  the  same 
way  as  described  previously  for  the  wing.  The  constant  is 
zero  below  the  forward  speed  input  in  XEL09  and  is  equal  to 
XEL08  above  the  forward  speed  input  in  XEL10.  Between  these 
two  speeds,  the  induced  velocity  effect  increases  linearly. 

For  helicopters  with  a  single  main  rotor,  the  speeds  usually 
used  are  zero  and  one  knot  so  that  the  elevator  is  always  in 
the  rotor  downwash.  For  a  helicopter  like  the  XV-3,  or  other 
configurations  which  use  laterally-disposed  rotors,  the  velo¬ 
cities  used  depend  upon  wake  skew  angles  and  the  location  of 
the  elevator.  The  multiplying  factor  usually  used  is  two, 
since  the  induced  velocity  which  is  computed  is  considered  at 
the  rotor  disk. 
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3.  Vertical  Stabilizer 

The  primary  geometric  and  aerodynamic  data  used  to  describe 
the  vertical  stabilizer  characteristics  are  similar  to  those 
described  previously  for  the  wing.  Two  aerodynamic  inter¬ 
ference  effects  are  simulated:  tail  rotor  induced  velocity  and 
the  rate  of  change  of  fin  angle  of  attack  with  sideslip 
(do/d p  in  conventional  terminology). 

The  tail  rotor  induced  velocity  is  multiplied  by  XFN06  and 
added  to  the  sideward  velocity  component  of  the  free  stream 
velocity  vector  acting  on  the  vertical  stabilizer.  The  tail 
rotor  mast  is  assumed  to  be  perpendicular  to  the  X-Z  plane 
of  the  aircraft.  If  the  vertical  stabilizer  is  to  the  right 
of  a  tail  rotor  which  is  providing  a  thrust  to  the  right 
(’’above"  the  rotor),  then  XFN06  is  usually  about  0.5.  If  the 
vertical  stabilizer  is  to  the  left,  then  XFN06  is  usually 
about  1.5,  depending  upon  proximity  of  the  fin  and  tail  rotor. 

The  free  stream  velocity  vector  is  modified  further  by  the 
sidewash  coefficient.  The  lateral  velocity  at  the  fin  is 
given  by: 


V  =  V  (1  -  XFN07) 
yfin  ybody 

If  no  sidewash  effect  is  desired,  then  XFN07  is  set  to  zero. 
Most  single  rotor  helicopters  have  a  value  of  about  0.2  for 
XFN07 . 

E.  JETS 

Auxiliary  propulsion  may  be  added  in  any  specified  amount  by 
means  of  the  input  data  of  the  jet  group.  Either  one  or  two 
jets  may  be  used. 

The  location  of  the  right  jet  is  specified  as  well  as  the 
pitch  and  yaw  angle  with  respect  to  the  body  axis  system. 

The  left  jet  is  assumed  to  be  at  the  same  station  and  water 
line  and  pitch  angle.  The  buttline  and  yaw  angle  of  the  left 
jet  are  assumed  to  be  the  negative  of  the  right  jet.  If  only 
one  jet  is  used,  it  is  assumed  to  be  the  "right”  jet  and  a 
negative  buttline  must  be  used  if  it  is  desired  to  locate  the 
jet  on  the  left  side  of  the  fuselage.  The  magnitude  of  jet 
thrust  does  not  vary  with  speed  in  the  current  program.  If 
no  jets  are  used,  the  number  of  jets  (XJET01)  is  set  to  zero. 

F .  CONTROL  SYSTEM 

The  synthesis  of  the  control  system  of  a  specified  aircraft 
configuration  is  complicated  because  of  the  generality  of  the 
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mathematical  model.  The  reader  is  referred  to  Volume  II  for 
a  detailed  discussion  of  synthesis  of  the  control  system. 

The  usual  control  system  linkages  used  for  three  typical 
configurations  are  indicated  in  Table  IX  by  x’s.  The  con¬ 
trolled  elements  are: 

Incidence  of  aerodynamic  surfaces: 

Horizontal  stabilizer 
Vertical  stabilizer 
Right  wing  panel 
Left  wing  panel 

Main  (forward,  right)  rotor  pitch  controls: 

Collective 
Longitudinal  cyclic 
Lateral  cyclic 

Tail  (aft,  left)  rotor  pitch  controls: 

Collective 
Longitudinal  cyclic 
Lateral  cyclic 

Auxiliary  thrust: 

Magnitude 

Direction 

These  elements  are  controlled  by  the  pilot  or  by  automatic 
systems.  The  pilot  controls  used  are: 

Collective  pitch  lever 

Longitudinal  cyclic  pitch  stick 

-  Lateral  cyclic  pitch  stick 

Pedals 

Mast  conversion  angle  lever 

The  automatic  systems  simulated  include  a  rotor  speed  governor, 
a  bobweight,  a  flat  tracker  mechanism,  and  focused  pylon 
geometry.  The  flat  tracker  equations  provide  cyclic  inputs 
which  act  in  such  a  way  as  to  minimize  rotor  flapping.  The 
focused  pylon  equations  change  cyclic  pitch  proportional  to 
the  rotor  in-plane  (H  and  Y)  forces. 
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G.  FLIGHT  PATH  EQUATIONS  OF  MOTION 

The  derivation  of  the  equations  of  motion  of  the  aircraft  is 
given  in  matrix  and  vector  form  in  Reference  4.  No  small 
angle  assumptions  are  made  in  the  analysis.  Euler  angles^^) 
are  used  to  orient  the  aircraft  with  respect  to  the  fixed 
reference  system.  Euler  angle  matrix  transformations  are  used 
to  transform  positions,  velocities,  and  accelerations  from  one 
axis  system  to  another  as  required. 

The  center  of  gravity  ,  gross  weight,  mass  moments,  and  products 
of  inertia  are  assumed  in  this  analysis  to  include  rotor  mast 
and  rotor  blade  weights.  The  change  in  center  of  gravity 
with  mast  conversion  angle  and  with  blade  fold  angle  is 
accurately  represented. 

The  rotors  are  assumed  to  be  coupled  together  at  a  fixed 
gearing  ratio.  This  enables  the  drive  system  to  be  approxi¬ 
mated  by  one  equation  involving  drive,  system  inertia,  supplied 
torque,  and  required  torque.  The  rotor  moments  about  the 
flapping  hinge  are  resolved  into  sine  and  cosine  components 
of  azimuth  position.  These  components  then  determine  the 
resulting  angular  accelerations  in  pitch  and  roll  of  the  rotor 
disk. 

A  summary  of  the  equations  is  given  in  Table  X. 
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SECTION  IV 


DETERMINATION  OF  THE  TRIM  CONDITION 

This  section  of  the  report  describes  some  of  the  pertinent 
details  of  the  method  used  to  compute  the  specified  trim  flight 
condition.  The  trim  solution  of  the  equations  must  be  deter¬ 
mined  before  any  of  the  other  capabilities  of  the  program  can 
be  used.  The  trim  solution  becomes  increasingly  difficult  as 
higher  flight  speeds  and/or  load  conditions  are  attempted. 

The  reason  for  this  is  that  the  terms  in  the  equations  of 
motion  become  more  nonlinear  as  the  effects  of  stall,  compres¬ 
sibility,  and  reverse  flow  become  more  predominant. 

A.  MATHEMATICAL  TECHNIQUE 

There  are  ten  equations  which  must  be  satisfied  by  the  trim 
subroutine.  The  six  equations  of  aircraft  forces  and  moments 
(X,Y,Z,L,M,  and  N)  are  presented  in  Table  X.  The  rotor  mo¬ 
ments  are  balanced  Ly  summing  the  contributions  of  all  blades 
and  then  setting  sine  and  cosine  components  of  the  following 
equation  to  zero  for  each  rotor.  The  equation  below  is 
obtained  by  rearrangement  of  the  blade  flapping  equation  on 
page  28.  The  flapping  positions  and  velocities  which  satisfy 
this  equation  then  determine  the  rotor  forces  and  moments. 

Pj  +  22  Ccosp  +  Tj)  sinp  =  0 


The  trim  solution  is  obtained  by  changing  the  values  of  the 
ten  variables  whose  values  are  unknown  and  which  affect  the 
above  equations.  These  variables  are  the  four  pilot  controls, 
the  longitudinal  and  lateral  flapping  angles  on  both  rotors, 
the  Euler  pitch  angle,  and  the  tenth  variable  is  either  the 
Euler  yaw  angle  or  roll  angle.  The  initial  estimates  of  these 
variables  which  are  part  of  the  input  data  (the  Flight  Con¬ 
stants  Group,  XFC)  must  be  reasonably  close  to  the  final  values 
in  order  for  the  trim  solution  procedure  to  work.  For  normal, 
level -flight  trim  conditions,  it  is  recommended  to  specify  all 
pilot  controls  at  50%,  all  flapping  angles  and  Euler  yaw  or 
roll  angle  at  zero,  and  use  some  reasonable  pitch  attitude 
depending  upon  the  configuration  and  speed.  The  values  of  the 
ten  variables  which  are  specified  in  the  input  data  are  then 
used  in  the  mathematical  model  to  evaluate  the  aerodynamic 
forces  and  moments.  If  the  error  between  the  desired  and 
computed  forces  and  moments  exceeds  the  allowable  values 
specified  in  the  Allowable  Error  Group,  XER,  then  changes 
must  be  made  to  the  values  of  the  variables  so  as  to  reduce 
the  error. 
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In  order  to  compute  the  amount  of  the  changes,  the  effect  of 
each  variable  on  the  ten  equations  is  evaluated  separately* 
This  is  done  by  incrementing  each  variable  by  a  small  amount 
(which  is  determined  by  the  data  supplied  in  the  Iteration 
Limits  Group,  XIT)  and  computing  the  resulting  changes  in  the 
forces  and  moments.  These  changes  are  then  divided  by  the 
amount  of  the  increment  of  the  variable  and  the  quotient 
represents  the  partial  derivative.  Mathematically,  this  pro¬ 
cedure  is  expressed  as  follows.  Let  the  equations  (i)  be 
represented  by  F  and  the  variables  (j)  by  x,  thus 


=  fCx^)  where  i,j  =  1  to  10 

For  the  increment  to  the  jth  variable,  the  ith  equation  can  be 
represented  as  the  sum  of  the  un- incremented  F^  and  the  change 
in  F^. 


Fi  +  AFi  =  f^xl  x2****xj  +  Axj'-*‘-xio) 
The  partial  derivative  is  then 


dF±  AFt 
= 

There  are  100  partial  derivatives  evaluated  and  they  are 
arranged  in  the  form  of  a  10x10  matrix.  When  the  errors  in 
forces  and  moments  are  determined,  these  ten  simultaneous 
equations  are  solved  for  the  changes  in  the  variables  which 
are  required.  Checks  are  then  made  to  make  sure  that  the 
magnitudes  of  the  computed  changes  are  not  too  large.  The 
new  values  of  the  variables  are  calculated  by  adding  the 
computed  changes  to  the  values  used  in  the  last  iteration. 

The  forces  and  moments  are  recomputed  and  a  new  partial  deriva¬ 
tive  matrix  is  determined  about  the  new  trim  point.  If  the 
equations  were  perfectly  linear,  there  would  be  no  change  in 
the  partial  derivative  matrix  between  iterations  and  only  one 
iteration  would  be  required.  This  is  rarely  the  case,  however, 
because  of  the  nonlinear  character  of  the  rotor  equations  and 
aerodynamic  surfaces  which  are  stalled.  An  option  is  available 
to  the  user  which  reduces  the  computer  time  required  to  trim 
in  linear  flight  regimes  or  when  fairly  small  changes  are 
being  made  by  the  sweep  feature  of  the  program.  This  option 
is  specified  by  a  non-zero  input  in  XIT3  which  then  causes  the 
program  to  recompute  the  partial  derivative  matrix  on  every 
fifth  iteration. 
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B.  CONVERGENCE 


The  above  described  procedure  usually  converges  to  a  trim 
solution  within  5  to  20  iterations.  Usually,  the  maximum 
number  of  iterations  allowed,  XIT1,  is  40.  For  some  config¬ 
urations  and  flight  conditions,  it  may  be  difficult  for  the 
above  technique  to  obtain  the  trim  solution  because  of  the 
nonlinearity  of  the  variables.  Mathematical  science  in  the 
area  of  solution  of  nonlinear,  simultaneous  equations  has  not 
advanced  to  the  point  where  solutions  can  always  be  found.  In 
fact,  most  techniques  will  fail  to  converge  unless  the  initial 
value  is  quite  close  to  the  final  solution. 

There  are  many  reasons  why  the  computer  program  may  fail  to 
converge  to  a  trim  solution.  Because  of  the  complexity  of  the 
program  as  well  as  the  mathematical  difficulties  described 
above,  it  is  not  possible  to  determine  a  list  of  causes  and 
cures  which  will  always  enable  a  trim  solution  to  be  found. 

The  first  check  to  make  if  trim  is  not  obtained  is  to  review 
all  of  the  input  values  for  keypunch  or  data  errors.  Described 
below  are  some  of  the  more  common  sources  of  trim  difficulties 
and  the  methods  for  getting  a  trim  solution. 

1 •  Iterative  Repetition 

Some  cases  will  use  all  of  the  allowed  number  of  iterations  and 
still  fail  to  trim.  For  these  cases,  the  last  sets  of  iteration 
data  should  be  compared  to  each  other  to  see  if  the  same  trim 
point  is  being  computed  every  other  (or  every  3rd,  4th,  etc.) 
time.  To  determine  whether  or  not  this  is  true,  the  values  of 
the  solution  variables  (VAR(I),  I  =  1  to  10)  are  compared  to 
see  if  they  are  all  equal.  The  most  common  causes  of  this 
problem  are  induced  velocity  or  airfoil  stall.  The  values  of 
X  — ,  Y-,  and  Z-  force  and  L-,  M-,  and  N-  moment  should  be  exam¬ 
ined  at  each  iteration  to  determine  the  magnitude  of  the  force 
changes  between  trim  conditions.  The  value  of  induced  velocity 
should  be  examined  for  both  rotors  to  see  if  changes  are  occur¬ 
ring  between  iterations  which  are  large  enough  to  induce  force 
changes  on  adjacent  aerodynamic  surfaces  that  are  a  significant 
fraction  of  the  total  force  change  between  iterations.  For 
example,  the  change  of  wing  load  caused  by  a  change  in  main 
rotor  induced  velocity  or  the  interaction  between  tail  rotor 
induced  velocity  and  fin  load.  If  this  occurs,  then  the  input 
data  specifying  the  appropriate  rotor  wake  effect  may  be  in 
error  or,  if  these  data  are  correct,  the  damping  term  on  in¬ 
duced  velocity,  XIT7 ,  should  be  decreased.  Normal  values  of 
XIT7  vary  from  0.5  to  1.0  and  may  have  to  be  decreased  to  0.1 
to  0.5  if  this  problem  occurs.  Sometimes  a  change  of  about  5 
degrees  of  Euler  angle  about  the  appropriate  axis  will  help 
get  out  of  the  induced  velocity  region  where  iterative  repeti¬ 
tion  occurs.  This  may  be  tried  if  the  above  changes  are  un¬ 
successful  . 
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Airfoil  stall  can  cause  iterative  repetition  if  one  trim  condi¬ 
tion  results  in  unstalled  flow  on  a  fixed  aerodynamic  surface 
and  the  following  trim  condition  results  in  stalled  flow.  The 
reason  for  this  is  the  change  in  sign  of  the  partial  derivative 
of  force  with  respect  to  angle  of  attack  which  is  caused  by  the 
change  in  the  slope  of  the  lift  curve  at  stall.  If  this  occurs, 
the  input  Euler  attitude  angles  can  be  changed  about  the  appro¬ 
priate  axes  to  unstall  the  surface.  If  a  high  CL  condition  is 
inherent  to  the  flight  condition  being  evaluated,  the  maximum 
Cli  YXX3,  may  be  input  at  a  larger  than  normal  value  in  order 
to  trim.  The  resulting  CL  can  be  compared  to  a  realistic  value 
to  determine  the  validity  of  the  results.  Another  source  of 
useful  information  about  the  above  problems  is  the  partial 
derivative  matrix  which  should  be  computed  and  printed  for 
every  iteration  if  trouble  is  encountered  in  trimming  the  air¬ 
craft  (XIT3  =  0) .  The  values  in  the  matrix  should  be  compared 
from  iteration  to  iteration  to  determine  significant  changes 
in  magnitude  and  signs  of  the  derivatives.  Suspect  derivatives 
can  then  be  further  traced  back  through  use  of  the  force  and 
moment  summaries  to  determine  likely  causes  of  the  observed 
changes.  This  will  often  indicate  what  changes  are  required  to 
obtain  a  trim  solution. 

2.  Excessive  Control  Travel 


For  some  cases,  the  iterative  technique  will  result  in  one  of 
the  controls  (VAR(I),  I  =  1  to  4)  exceeding  100%  or  less  than 
0%  travel.  Trim  may  or  may  not  be  obtained  in  this  instance. 

If  trim  is  obtained,  then  an  examination  of  the  trim  solution 
will  usually  indicate  which  aircraft  parameters  must  be  altered 
to  bring  the  control  travel  into  the  normal  range. 

If  trim  is  not  obtained,  then  the  force  and  moment  summary  data 
at  each  iteration  should  be  examined  to  determine  which  ele¬ 
ments  are  causing  excessive  forces  and  moments.  The  cause  of 
the  increase  in  the  force  or  moment  of  the  element  should  be 
isolated  if  possible.  The  partial  derivative  matrix  should 
then  be  examined  to  determine  which  are  the  primary  derivatives 
(indicated  by  largest  numerical  value)  relating  variables  to 
the  excessive  force  or  moment.  These  variables  should  then  be 
changed  in  the  appropriate  direction,  indicated  by  the  sign  of 
the  partial  derivative,  and  then  rerun  the  case. 

3.  First  Iteration  Technique 

If  the  methods  described  above  should  fail  to  attain  a  trim 
condition,  there  is  a  method  for  enabling  the  user  to  better 
understand  the  dependency  of  the  forces  and  moments  on  the 
initial  values  which  are  used  in  the  input  data. 
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The  input  values  of  allowable  force  and  moment  errors,  XER1  to 
XER7 ,  are  set  very  large  (9,000,000)  so  that  the  computer  recog 
nizes  a  "trim"  solution  no  matter  what  the  input  data  values 
are.  All  data  are  then  printed  out  on  both  the  iteration  and 
trim  page.  Secondly,  the  sweep  feature  is  used  to  enable  the 
"trim"  data  to  be  computed  for  a  suitable  range  of  the  input 
parameters  which  are  specified  by  the  sweep  cards.  The  result¬ 
ing  variations  in  forces  and  moments  with  changes  in  the  input 
variables  are  noted  and  a  better  approximation  of  the  values 
of  the  variables  for  trim  can  be  made.  By  using  values  which 
are  closer  to  the  trim  condition,  convergence  of  the  computer 
program  to  the  trim  solution  is  more  likely.  A  major  advantage 
of  this  technique  is  the  extremely  short  amount  of  time  needed 
to  compute  data  for  one  iteration. 
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SECTION  V 


DISCUSSION  OF  TYPICAL  APPLICATIONS 


The  purpose  of  this  section  is  to  describe  some  of  the  more 
common  applications  of  the  computer  program  and  to  assist  the 
user  in  gaining  an  appreciation  of  the  effects  of  changes  in 
the  values  of  some  of  the  input  data  on  the  computed  flight 
characteristics.  One  of  the  major  obstacles  to  the  use  of 
the  computer  program  is  the  time  and  effort  required  to  deter¬ 
mine  the  values  for  over  400  items  of  input  data.  Once  this 
is  done,  however,  a  powerful  analytical  tool  is  available 
which  has  extensive  application  during  preliminary  design, 
flight  testing,  and  for  investigating  the  effects  of  aircraft 
geometry  changes,  weapon  recoil,  external  stores,  and  a  vari¬ 
ety  of  other  parameters  on  the  flight  characteristics. 

It  cannot  be  overemphasized  that  the  qual¬ 
ity  of  the  results  will  never  exceed  the 
quality  of  the  input  data. 

The  mathematical  model  in  the  computer  program  is  very  real¬ 
istic.  If  the  input  data  define  an  unattainable  flight  condi¬ 
tion  because  of  physical  laws  or  power  limitations,  then 
either  no  trim  condition  will  be  found  or  some  aspect  of  the 
output  data  will  not  be  realistic.  The  user  must  recognize 
this  condition  and  correct  the  input  data.  Some  of  the  input 
data  require  experience  and  judgment  in  the  determination  of 
values  to  use.  These  date  are  discussed  below  to  assist  the 
user  in  acquiring  a  working  knowledge  of  the  program. 

A.  PERFORMANCE  CORRELATION 

This  is  generally  the  first  topic  of  interest  after  the  initial 
values  for  input  data  have  been  determined  from  the  geometric 
and  physical  characteristics  of  the  aircraft.  The  extent  of 
the  correlation  will  be  determined  by  the  purpose  of  the  sub¬ 
sequent  investigations  as  well  as  the  amount  and  reliability 
of  data  to  which  comparison  will  be  made.  Data  measured  in 
flight  test  will  probably  be  best  if  they  have  been  carefully 
measured  and  reduced.  In  the  design  phase  of  an  aircraft, 
predictions  made  by  the  Aerodynamics  or  Performance  Group  are 
generally  used. 

1 .  Power  Required 

The  total  power  required  is  computed  for  all  trim  flight  con¬ 
ditions.  The  parameter  sweep  feature  may  be  used  to  change 
gross  weight,  speed,  and  drag  area,  while  using  only  one  basic 
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set  of  input  data.  Drag  area  changes  must  often  be  considered 
because  the  configuration  of  external  stores  will  change  the 
effective  drag  area  of  the  aircraft. 

The  major  contributors  to  rotorcraft  power  required  are: 
flat-plate  drag  area 
rotor  induced  power 
rotor  profile  power 

There  are  several  input  data  that  have  a  first-order  effect  on 
the  magnitude  of  each  of  the  above  terms.  The  magnitude  of 
these  terms  varies  with  speed  and  weight  in  very  different 
manners.  These  relationships  will  be  reviewed  to  assist  the 
user  in  determining  which  power  contributor  should  be  modified 
to  improve  correlation. 

In  the  computer  program,  flat-plate  drag  area  is  divided  into 
the  following  primary  components. 


Component 

Input  Data  Specified 

fuselage 

XFS22  through  XFS25 

wing 

XWG1  and  YWG12  to  YWG14 

elevator 

XEL1  and  YEL12  to  YEL14 

fin 

XFN1  and  YFN12  to  YFN14 

jets 

XJET2  and  XJET3 

rotor  hub 

XMR21  and  XMR25 

For  initial  performance  work,  the  drag  contribution  of  the 
above  components  is  not  as  important  as  the  total  flat-plate 
drag.  Reasonable  estimates  of  drag  coefficients  and  areas 
are  adequate  at  this  point.  The  effect  of  flat  plate  drag 
on  power  required  is  a  velocity-cubed  function.  It  is  the 
product  of  the  drag  force  and  the  velocity  cubed. 

Pf  =  0 . 5pfV3 

Knowledge  of  this  characteristic  relationship  is  useful  when 
making  small  corrections  to  the  input  data,*  above,  in  order  to 
reduce  power  discrepancies  which  vary  as  V3. 
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Rotor  induced  power  is  determined  internally  and  is  not  sub¬ 
ject  to  modification  by  altering  input  data  values.  Induced 
power  can  be  approximated  by  the  product  of  the  rotor  thrust, 
usually  equal  to  the  gross  weight,  and  an  induced  velocity 
which  is  based  upon  momentum  theory.  This  induced  velocity 
is  calculated  by  solving  the  following  equation  for  VD  in  an 
iterative  manner. 


v  -  =  GW 

p  GW  2pA  /Vp2  +  V2 

where  represents  the  flat  plate  drag  power  discussed 
previously. 

Once  V  is  determined,  it  is  used  in  either  side  of  the  above 
equation  which  is  equal  to  the  induced  velocity,  i.e., 

pf 

vi  =  VP  -  5w 

In  hovering  flight,  it  is  apparent  that  the  above  equation 
for  reduces  to 


V. 

1 


2  2 

and  in  forward  flight  above  about  60  knots,  V  »Vp 


Vi 


^  GW 
-  Tpffl 


so 


These  approximate  expressions  are  often  adequate  for  correla¬ 
tion  purposes.  The  actual  induced  velocity  equations  used  are 
much  more  complex  and  account  for  spanwise  variation  and 
ground  effect. 


Rotor  profile  power  is  nearly  constant  but  varies  somewhat  with 
thrust,  or  gross  weight,  as  well  as  speed.  The  variation  with 
speed,  which  is  caused  by  flow  asymmetry  in  forward  flight,  is 
inherent  in  the  blade  element  representation  of  the  rotor  and 
is  not  subject  to  user  modification  by  means  of  input  data. 

The  variation  with  thrust  is  also  built  into  the  rotor  repre¬ 
sentation  but  may  be  modified  by  the  user  through  the  input 
data  for  the  airfoil  characteristics  subroutine.  The  approxi¬ 
mate  effect  on  power  of  varying  input  data  in  the  rotor  aero¬ 
dynamic  group  will  be  discussed  below.  The  exact  drag  equation 
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used  is  described  in  Section  IIIB  and  includes  stall  and 
compressibility  effects. 


Most  of  the  rotor  operates  in  the  unstalled ,  subsonic  flow 
regime.  The  drag  coefficient  equation  for  a  main  rotor  blade 
element  operating  in  this  regime  may  be  approximated  by 

CD  =  YMR12  +  a(YMR13  +  aYMR14) 

where  a  is  the  local  section  angle  of  attack.  For  purposes 
of  overall  correlation,  an  average  value  of  a  is  needed. 

This  value  may  be  arrived  at  through  the  expression  for 

and  the  input  YMR17.  Since,  for  symmetrical  airfoil  sections, 

YMR13  should  be  zero,  the  following  may  then  be  used: 


or,  in  dimensional 
units , 


/  7ct  V 

\amcsj) 


YMR12  ♦ 

terms,  with  YMR17  and  YMR14  in  consistent 

2 


YMR12  ♦  YMR14 


( 


7  GW 


p(YMR17)bcR(SR) 


: 


where  b  is  the  number  of  blades.  The  profile  power  may  then 
be  approximated  by: 


P0  =  CD 


(f)pA(2R)3  1  +  4.6  (^) 


The  effectiveness  of  the  inputs  YMR12  and  YMR14  on  power  re¬ 
quired  is  apparent  after  substitution  into  the  preceding  two 
equations.  The  input  values  can  then  be  adjusted  within 
reasonable  limits  to  improve  power  correlation  with  variation 
in  speed  and  gross  weight . 


Power  divergence  is  sometimes  apparent  in  test  data  above  a 
certain  forward  speed.  This  is  caused  by  compressibility 
effects  on  drag.  These  effects  also  appear  as  different 
Cp-Cx  curves  with  different  rotor  speeds.  The  input  variable 
which  affects  the  magnitude  of  compressibility  effects  is 
YMR1,  the  critical  Mach  number  of  the  airfoil  used.  YTR1 
affects  tail  rotor  power  in  a  similar  manner.  Typical  values 
for  YM^l  vary  from  .78  to  .88  depending  primarily  on  thickness. 

2.  W^ng-Rotor  Lift  Sharing 

One  of  two  problems  are  generally  considered  in  this  regard. 
Either  correlation  is  desired  between  measured  and  computed 
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data,  or  geometric  and  control  system  parameters  are  being 
optimized  during  the  design  phase  to  achieve  a  specified 
lift  sharing  ratio. 

In  the  former  case  of  correlation  with  measured  data,  only 
one  input  will  directly  alter  lift  sharing  significantly. 

This  input  is  XWG8  which  multiplies  the  induced  velocity  of 
the  main  rotor  and  adds  the  resulting  velocity  vector,  which 
acts  along  the  mast,  to  the  free  stream  velocity  vector  to 
obtain  the  wing  angle  of  attack.  Typical  values  of  XWG8  are 
1.0  to  2.0,  depending  upon  wing  location.  Another  significant 
factor  is  the  variation  of  fuselage  pitch  attitude  with  for¬ 
ward  speed,  gross  weight,  center  of  gravity  location,  and 
external  stores  configuration.  All  inputs  which  have  an  effect 
on  pitching  moment  equilibrium  will  affect  the  fuselage  pitch 
attitude.  Listing  all  of  these  inputs  would  take  too  long  and 
be  too  dependent  on  the  configuration  to  be  valid.  Also,  many 
of  the  terms  which  affect  pitching  moment  are  determined  by 
other  considerations.  Some  of  these  terms  are  hub  spring 
restraint,  elevator  incidence  and  location  and  gearing,  and 
wing  incidence  and  location.  It  is  desired  to  find  input 
values  that  can  be  altered  within  reasonable  limits  to  improve 
correlation  and  values  which  do  not  have  strong  impact  in  other 
areas  of  the  design.  These  inputs  are  usually  the  fuselage 
pitching  moment  coefficients,  XFS15  and  XFS16.  Fuselage  lift 
and  drag  characteristics  and  the  location  of  the  aerodynamic 
center  may  be  altered  somewhat,  but  these  changes  will  sig¬ 
nificantly  affect  cyclic  control  stick  margins  and  gradients 
as  discussed  in  the  subsequent  section. 

If  the  desired  lift  sharing  characteristics  are  known  for  a 
specified  design,  the  user  has  freedom  to  determine  geometric 
and  control  system  parameters  in  order  to  meet  the  desired 
lift  schedule.  As  these  parameters  are  changed,  the  effects 
on  stability,  control  margins,  control  gradients,  and  power 
must  also  be  considered.  A  combination  of  parameters  that 
satisfies  the  lift  schedule  at  the  expense  of  violating  one 
of  the  other  areas  mentioned  is  of  little  practical  use. 

If  the  wing  has  no  movable  surfaces  or  is  not  wholly  mov¬ 
able,  then  the  control  system  parameters  will  not  have  a 
first-order  effect  on  lift  sharing.  If  the  wing  incidence  is 
variable,  the  effect  of  different  combinations  of  changing 
incidence  with  collective  pitch  and/or  longitudinal  cyclic 
pitch  may  be  determined  by  changing  the  appropriate  inputs 
in  the  Controls  Group,  the  XCON  inputs.  Geometric  parameters 
of  the  wing  that  are  usually  altered  to  change  the  lift  sched¬ 
ule  are:  area,  XWG1;  horizontal  location,  XWG2;  incidence, 
XWG5;  and  to  some  extent,  aspect  ratio,  YWG18.  The 
first  two  of  these  parameters  have  a  strong  influence  on 
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control  margins,  gradients,  and  aircraft  stability  so  these 
effects  must  be  examined  concurrently  with  the  lift -sharing 
schedule . 

B.  STABILITY  AND  CONTROL  ANALYSIS 

Control  positions  are  determined  for  all  trim  flight  condi¬ 
tions.  The  control  power  derivatives  are  also  determined 
which  relate  the  change  of  pitching,  rolling,  and  yawing 
moment  to  control  deflection.  The  stability  analysis  and 
vector  analysis  sections  may  also  be  used  to  determine  the 
effects  of  parameter  changes  on  the  flying  qualities  of 
interest . 

The  computer  program  is  usually  used  either  to  correlate 
measured  data  or  to  determine  geometric  and  control  system 
parameters  of  new  aircraft  designs.  The  user  generally  has 
more  latitude  in  altering  values  of  input  data  during  design 
studies  than  for  correlation  studies  of  actual  aircraft. 
Although  this  gives  greater  flexibility,  it  also  increases 
the  number  of  variables  being  considered,  so  care  must  be 
taken  in  selecting  which  variables  to  alter  and  in  what 
combinations . 

1 .  Control  Margins  and  Gradients 

After  performance  correlation  is  obtained,  the  control  posi¬ 
tion  and  gradient  variation  with  changes  of  airspeed,  gross 
weight,  center  of  gravity  location,  and  external  store  con¬ 
figuration  are  then  examined.  If  the  correlation  is  not  ade¬ 
quate,  some  of  the  input  data  may  be  adjusted  within  reason¬ 
able  limits  to  improve  correlation. 

The  amount  of  longitudinal  cyclic  control  which  is  required 
to  trim  pitching  moments  will  vary  with  forces  and  moments 
acting  on  the  fuselage.  The  variables  which  affect  this  are 
the  same  as  those  which  affect  pitch  attitude.  In  addition 
to  these,  the  apportionment  of  drag  area  between  the  fuselage 
and  the  hub  may  be  altered  somewhat  because  of  the  uncer¬ 
tainty  associated  with  determining  an  accurate,  full-scale  hub 
drag  value.  Re -proportioning  the  drag  areas  will  change  the 
pitching  moment  acting  oi.  the  fuselage  and  result  in  a  change 
in  control  stick  position  required  to  trim.  The  control 
power  data  are  useful  in  determining  how  much  of  a  change  is 
required  in  pitching  moment  to  move  the  control  a  specified 
amount.  Since  the  total  drag  area  does  not  change,  the  power 
correlation  will  not  be  affected  significantly. 

Lateral  cyclic  and  pedal  control  required  to  trim  will  vary 
with  forces  and  moments  acting  on  the  fuselage  and  fin.  In 
the  absence  of  wind-tunnel  data,  several  variables  may  be 
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varied  within  reasonable  limits  because  of  the  difficulty  of 
measuring  or  estimating  their  effects  precisely.  These 
variables  are  : 

effective  fin  aspect  ratio,  YFN18 

-  effective  fin  area,  XFN1 

fin  aerodynamic  center,  XFN2  and  XFN4 

sidewash  coefficient,  XFN7 

Fuselage  yaw  moment ,  XFS17  and  XFS18 

Fuselage  side  force,  XFS26  to  XFS28 

Fuselage  aerodynamic  center,  XFS2  and  XFS4 

Although  data  obtained  in  the  wind  tunnel  are  better  than 
estimated  data,  exact  correlation  will  probably  not  be 
achieved  because  of  Reynolds  Number  effects  and  the  effect  of 
the  actual  rotor  flow  field  impingement  on  the  fuselage  and 
aerodynamic  surfaces.  Some  latitude  may  be  used  in  the  values 
assigned  to  the  input  variables  to  account  for  these  effects. 

2 .  Control  System  Kinematics 

A  wide  variety  of  control  system  linkages  may  be  synthesized 
by  input  data,  as  previously  described  in  Section  IIF  and 
Table  IX.  This  capability  is  especially  useful  when  design¬ 
ing  the  control  system  of  new  aircraft.  Elevator  synchroni¬ 
zation  schedules  with  collective  pitch,  and  longitudinal 
cyclic  pitch  and  mast  tilt  may  be  altered  by  the  inputs  XC0N11, 
XCON27 ,  XCON28,  XCON63  and  XEL5 .  The  corresponding  changes 
in  control  margins  and  gradients  are  then  computed  for  the 
specified  configuration  (weight,  eg,  external  store)  and 
flight  condition  (speed,  load  factor,  altitude). 

The  tilt-rotor  concept  uses  differential  collective  pitch  for 
roll  control  in  helicopter  mode  and  differential  longitudinal 
cyclic  pilch  for  yaw  control.  As  the  rotors  are  converted 
through  ninety  degrees,  these  swashplate  motions  result  in 
lateral  cyclic  control  producing  yaw  motion  and  pedal  control 
producing  a  roll  motion  unless  some  means  of  phasing  these 
controls  with  conversion  angle  are  used.  In  addition  to 
this  change  of  control  response,  the  rotor  cyclic  and  differen¬ 
tial  collective  controls  must  be  phased  out  so  that  only 
collective  pitch  is  controlled  by  the  pilot  (or  governor)  in 
airplane  mode.  In  order  to  investigate  how  these  controls 
should  be  phased,  the  computer  program  was  used  to  compute 
yaw  and  roll  moment  caused  by  pedal  deflection  with  different 
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amounts  of  differential  collective  pitch  being  actuated  by 
pedal  motion.  The  results  are  summarized  in  Figure  6.  It 
was  desired  to  determine  the  coupling  that  would  minimize  the 
roll  moment  caused  by  pedal  deflection  while  maintaining  a 
linearly  decreasing  value  of  yaw  control  power.  The  dashed 
lines  achieve  these  results  at  the  design  conversion  speed  if 
a  15°  mast  tilt  is  used  above  80  to  100  knots.  The  resulting 
control  kinematics  are  also  shown  in  Figure  6.  The  effects  of 
other  control  system  linkages  for  other  aircraft  or  helicopter 
configurations  may  be  easily  determined  in  a  similar  manner. 
Other  topics  of  interest  may  include  wing  flap/aileron  linkage, 
rudder-tail  rotor  lift-sharing,  elevator  synchronization, 
other  control  phasing  with  mast  conversion,  effects  of  gover¬ 
nors  and  bobweights  and  other  parameters. 

3.  Stability  Characteristics 

The  stability  derivatives  are  computed  by  the  increment 
technique  discussed  in  detail  in  Section  IIB.  The  contribu¬ 
tions  of  the  various  components  of  the  aircraft  to  any  sta¬ 
bility  derivative  may  be  determined  by  examination  of  the  out¬ 
put  data  as  shown  in  Table  I.  The  accuracy  of  the  derivatives 
may  be  rapidly  determined  from  the  data  available  and  by  using 
the  standard  equations  for  the  fuselage  and  fixed  aerodynamic 
surfaces.  The  input  parameters  that  alter  these  derivatives 
will  be  apparent  to  the  fixed-wing  aerodynamicist . 

The  rotor  force  and  moment  derivatives  with  respect  to  linear 
and  angular  rates  will  correlate  well  in  low-speed  flight  with 
derivatives  which  are  available  in  the  literature. ( 15 In 
high-speed  flight,  the  more  accurate  representation  of  the 
rotor  in  the  computer  program  will  compute  more  accurate  sta¬ 
bility  derivatives.  These  derivatives  will  include  the  non¬ 
linear  effects  of  stall,  reverse  and  radial  flow,  and  compres¬ 
sibility  which  are  not  subject  to  simple,  closed-form  analysis. 

The  rotor  derivatives  may  not  be  computed  in  the  stability 
section  if  the  relationship  between  the  increment  size,  XIT4, 
used  and  the  allowable  error  in  rotor  moment  balance,  XER6  and 
XER7,  is  not  properly  proportioned.  The  reason  for  this  is 
that  the  rotor  subroutine  will  only  compute  new  flapping  angles 
if  the  linear  or  angular  rates  cause  a  change  in  flapping 
moment  that  exceeds  the  allowable  error  in  flapping  moment. 

If  the  rate  changes  are  too  small  or  if  the  allowable  moments 
are  too  high,  then  no  change  in  flapping  will  occur  because  of 
rates,  and  the  rotor  force  and  moment  derivatives  will  be  zero. 
This  condition  is  illustrated  in  Figure  2  where  the  increment 
size,  XIT4,  is  below  a  value  of  0.04;  that  is,  the  change  in 
linear  velocities  is  0.4  feet  per  second  and  the  change  in 
angular  velocities  is  0.004  radians  per  second.  The  usual 
value  of  XER6  that  is  used  is  50  pound-feet,  which  was  the 
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value  used  in  this  140  knot  sample  case.  A  velocity  increment 
of  0.4  feet  per  second  perpendicular  to  the  140  knot  velocity 
vector  will  induce  an  angle  of  attack  or  side-slip  change  of 
only  .0017  radians.  This  small  a  change  is  insufficient  to 
cause  the  program  to  rebalance  the  flapping  moment.  It  should 
also  be  noted  that  for  very  low  increment  sizes  the  resulting 
small  differences  of  large  numbers  will  cause  a  wide  variance 
in  the  values  of  the  computed  derivatives.  A  good  rule  of 
thumb  for  the  increment  size  to  use  in  the  stability  analysis 
is  XIT4  =  0.002V,  where  V  is  in  feet  per  second.  For  the  sample 
case  of  Figure  2,  XIT4  should  be  0.5. 

If  the  program  is  being  used  for  correlation  with  measured 
data,  relatively  few  inputs  may  be  changed  to  improve  correla¬ 
tion.  Experience  indicates  that  correlation  between  measured 
and  computed  frequencies  and  damping  characteristics  is  good, 
and  the  input  data  seldom  needs  changing  if  the  initial  values 
were  carefully  determined.  Errors  in  input  data  can  often  be 
traced  back  by  means  of  the  stability  analysis  section  which 
provides  the  user  with  data  concerning  the  contribution  of 
each  element  to  a  stability  derivative.  Some  of  the  more 
important  input  data  in  this  regard  are  those  which  account 
for  aerodynamic  interference  effects.  These  data  are: 

Main  rotor  wake  on  wing 

Wing  wake  on  horizontal  stabilizer,  ri 

Main  rotor  wake  on  elevator 

Variation  of  XEL8  with  forward  speed 

Wing  wake  on  horizontal  stabilizer,  de/da 
Tail  rotor  wake  on  fin 

Fin  sidewash  caused  by  fuselage  and  wing 

The  program  may  be  used  to  compare  the  stability  characteristics 
of  a  new  aircraft  to  requirements  and  then  change  appropriate 
design  parameters  if  the  requirements  are  not  met.  In  this 
case,  best  approximations  of  the  above  interference  effects 
should  be  made  based  upon  the  specific  configuration.  Since 
the  rotor  induced  velocity  is  computed  at  the  disk,  the 
coefficient  to  use  to  simulate  an  airfoil  in  a  fully  developed 
rotor  wake  is  2.0,  but  values  of  0.5  to  1.0  are  more  represen¬ 
tative  of  actual  flow  conditions.  The  flow  angularity,  XEL  11 
and  XFN7 ,  at  the  empennage  which  is  caused  by  the  wing  may  be 
estimated  by  conventional  methods  (12,13,16)  if  wind  tunnel 
data  are  not  available. 

The  input  parameters  most  commonly  varied  to  study  control- 
fixed  stability  are  those  relating  to  the  size,  incidence. 


-  XWG8 

-  XWG9  l 

-  XWGIOf 

-  XEL8 

-  XEL9  \ 

-  XEL10J 

-  XEL11 

-  XFN6 

-  XFN7 
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aspect  ratio,  and  placement  of  the  fixed  aerodynamic  surfaces. 
The  sweep  technique  discussed  in  Section  IIA  enables  the 
effect  of  changes  of  these  parameters  on  stability  to  be 
easily  determined.  The  computed  frequency  and  damping 
characteristics  enable  direct  comparisons  with  specification 
requirements  to  be  made.  Since  trim  flight  data  may  be 
obtained  in  climbing  and  descending  flight,  the  effect  of 
power  on  control  position  may  be  directly  determined.  The 
neutral  point  may  be  determined  by  computing  stability 
derivatives  as  the  center  of  gravity,  XFS5,  is  moved  aft.  The 
eg  location  at  which  the  change  of  pitching  moment  with  angle 
of  attack,  dM/da,  is  zero  is  the  neutral  point.  The  maneuver 
margin  may  be  determined  in  a  similar  manner  by  examining  the 
change  in  the  longitudinal  cyclic  pitch  required  to  trim  as 
the  eg  is  moved.  The  eg  location  for  which  the  longitudinal 
cyclic  pitch  does  not  change  as  the  eg  is  moved  determines 
the  maneuver  point . 

Many  other  applications  of  the  program  to  stability  and  con¬ 
trol  investigations  are  possible.  The  examples  given  above 
are  intended  to  be  illustrative  and  not  exhaustive. 

G.  MANEUVERING  FLIGHT 

The  maneuver  computing  section  of  the  program  is  used  to 
graphically  illustrate  the  control  and  gust  response  character¬ 
istics  of  the  aircraft.  Many  flying  qualities  specifications 
are  expressed  in  terms  of  a  range  of  acceptable  deviations 
of  certain  parameters  during  the  performance  of  a  maneuver. 
Although  this  type  of  specification  is  intended  primarily  for 
application  in  flight  test,  the  maneuver  computing  section  of 
the  program  enables  the  maneuver  to  be  simulated  to  determine 
specification  compliance.  The  effect  of  design  parameters 
on  the  response  may  be  evaluated  for  new  designs,  or  the 
computed  and  measured  maneuver  time  histories  may  be  compared. 
If  the  measured  response  is  not  acceptable,  the  computer 
program  may  be  used  to  investigate  different  ways  of  obtaining 
acceptable  response. 

1.  Nonlinear  Effects 

Since  the  equations  of  the  maneuver  section  of  the  program 
are  fully  coupled  and  nonlinear,  response  to  simultaneous 
longitudinal  and  lateral  control  or  gusts  may  be  computed. 

The  vector  analysis  section  can  be  used  to  represent  as  a 
vector  quantity  any  parameter  that  can  be  plotted.  The 
magnitudes  and  phase  angles  of  the  vectors  can  then  be  compared 
to  the  mode  shapes  computed  in  the  stability  analysis  section. 
If  the  vector  relationships  and  mode  shapes  are  about  equal, 
as  they  usually  are,  then  the  nonlinear  effects  are  small 
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and  the  data  obtained  in  the  linear,  uncoupled  stability 
section  accurately  represents  the  stability  characteristics. 


The  correlation  between  these  methods  is  not  automatically 
assured.  Although  the  mathematical  model  is  the  same,  the 
techniques  used  to  compute  stability  are  quite  different. 

The  stability  section  uses  the  computed  stability  derivatives 
in  linear  perturbation  equations  and  then  computes  the  roots 
of  the  characteristic  equation.  The  maneuver  section  con¬ 
tinuously  extrapolates  to  the  next  maneuver  time  point  by 
applying  the  Runge-Kutta  technique  (described  in  Section  IIC) 
to  the  differential  equations  of  motion  of  the  aircraft,  a 
partial  list  of  which  is  in  Table  X. 

The  data  in  the  stability  section  are  usually  computed  in  less 
than  a  minute  on  the  IBM  360,  Model  65  computer.  Maneuver 
cases  take  approximately  one  minute  of  computing  time  per 
second  of  maneuver  time.  This  restricts  most  maneuver  invest i 
gations  to  short  period  phenomena  and  encourages  the  use  of 
the  stability  section  method,  since  most  flight  conditions 
are  characterized  by  linear  response. 

2.  Rotor  Flapping 

The  stability  section  computes  the  rotor  derivatives  with 
respect  to  flapping,  but  no  information  is  provided  concerning 
flapping  response  to  control  motion  or  gust  penetration.  The 
maneuver  section  is  well  suited  to  this  type  of  work  because 
the  high  frequency  and  damping  of  flapping  motion  require  a 
relatively  short  maneuver  time. 

The  amount  of  rotor  flapping  freedom  that  must  be  provided 
in  both  accelerated  and  unaccelerated  flight  may  be  determined 
from  the  data  computed  in  the  trim  section.  During  a  rotor 
design  study,  many  flight  conditions  are  examined  which 
determine  different  parameters  of  the  aircraft.  Once  the 
flight  conditions  which  exhibit  the  larger  flapping  angles 
have  been  determined,  the  maneuver  section  is  used  to  compute 
the  flapping  resulting  from  control  motions  and  gust  penetra¬ 
tion  in  the  direction  which  tends  to  increase  flapping. 

For  rotors  with  no  flapping  freedom,  the  mathematical  model 
provides  for  a  nonisentropic  hub  spring  of  specified  stiffness 
(XMR17  and  XMR18)  which  is  located  at  a  specified  radial 
location  (XMR16).  The  hub  moment  of  a  "rigid”  rotor  is  then 
indicated  by  the  product  of  the  computed  flapping  angle  and 
the  hub  restraint.  Changes  in  design  parameters  may  be 
made  in  the  input  data*,  and  the  effects  on  flappirg  will  be 
computed  for  the  design  flight  conditions.  In  the  event  that 
unacceptable  flapping  angles  or  hub  moments  are  computed, 
suitable  changes  must  be  made  in  the  design  to  reduce  flapping 
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Some  of  the  more  effective  parameters  which  reduce  flapping 
in  maneuvering  flight  are  those  which  reduce  rotor  blade 
loading  (increased  chord  or  radius)  or  reduce  the  amount  of 
thrust  change  required  from  the  rotor  in  order  to  attain  a 
specified  normal  acceleration.  Although  a  wing  is  effective 
in  altering  the  latter,  the  size  must  not  be  so  large  as  to 
unload  the  rotor  in  autorotation  (loss  of  rotor  speed  control) 
or  require  too  great  a  change  in  the  longitudinal  cyclic  pitch 
required  to  trim  pitching  moments  caused  by  power  changes. 

3.  Response  Characteristics 

Two  types  of  response  are  generally  of  interest:  control 
response  and  gust  response.  There  are  several  different 
shapes  of  gusts  that  may  be  simulated  in  the  computer  program 
as  described  previously  in  Section  IIC.  This  section  will 
discuss  some  of  the  techniques  used  when  investigating 
response  characteristics. 

Some  aircraft  exhibit  strong  coupling  between  pitching,  rolling, 
and  yawing  motion  because  of  the  inherent  configuration  of  the 
aircraft  or  because  of  the  flight  condition.  It  is  often 
desirable  to  study  the  effects  of  design  parameters  on  response 
about  one  axis  at  a  time  or  about  one  pair  of  axes  at  a  time. 
This  is  easily  simulated  in  the  computer  program  by  increasing 
the  inertia  about  the  axes  about  which  no  motion  is  desired. 

A  typical  value  for  these  inputs  (XFS8  to  XFS10)  to  lock  out 
motion  would  be  9,000,000  slug-feet^. 

When  the  program  is  used  to  compute  the  response  to  control 
inputs,  some  care  must  be  exercised  not  to  use  step  inputs. 

Not  only  are  step  inputs  impossible  to  attain  in  practice, 
they  also  cause  a  divergent  oscillation  to  occur  in  the  time- 
variant  solution  of  the  differential  equations  of  motion. 

This  divergence  would  not  occur  on  an  aircraft,  but  is  caused 
by  the  inability  of  the  Runge-Kutta  technique  to  accurately 
represent  step  changes  of  parameters.  When  pilots  are  requested 
to  make  "rapid"  or  "step"  control  inputs,  they  usually  use 
from  0.15  to  0.20  seconds  to  make  the  control  motion.  At 
least  this  amount  of  time  should  be  used  to  make  control 
position  changes  in  the  maneuver  section. 

If  a  step  input  must  be  made,  it  can  be  effectively  simulated 
by  using  a  very  small  control  input  duration,  on  the  order 
of  0.001  seconds,  and  a  time  increment  ( At  =  0.0001  seconds) 
which  results  in  at  least  ten  control  increments  to  accomplish 
the  input.  This  method  is  as  effective  as  a  step  input,  and 
avoids  the  divergence  resulting  from  a  large  step  change  in 
the  mathematical  model. 
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4.  Conversion  Maneuvers 

The  user  must  specify  the  position  of  all  pilot  controls  during 
the  computation  of  maneuver  response  characteristics.  Before 
a  conversion  maneuver  can  be  computed  then,  the  required  con¬ 
trol  positions  must  be  known  in  advance.  This  may  be  done 
by  computing  trim  flight  conditions  for  partial  conversion 
configurations  and  then  using  the  control  position  versus 
conversion  relationship  to  determine  the  required  time 
dependent  control  positions  during  the  conversion  maneuver. 

For  example,  consider  the  TFTA  aircraft  operating  in  heli¬ 
copter  mode  with  rotor  masts  vertical.  The  parameter  sweep 
feature  can  be  used  to  determihe  the  trim  control  positions 
at  15,  30,  45,  60,  75,  and  90  degrees  of  rotor  mast  conversion 
angle  at  the  same  airspeed.  After  these  data  have  been  com¬ 
puted,  they  can  be  plotted  as  a  function  of  time  once  the 
user  specifies  the  conversion  angle  versus  time  relationship. 
The  resulting  time  histories  of  control  motion  can  then  be 
approximated  by  a  series  of  control  rates  and  time  durations 
and  then  input  to  the  maneuver  section  on  appropriate  "J" 
cards  CVol.  II,  p.  60). 

This  same  technique  can  be  used  during  rotor  stopping  by 
computing  trim  conditions  at  different  rotor  speeds  as  well 
as  rotor  blade  fold  angles .  Care  must  be  used  to  insure 
that  all  pilot  controls  have  effective  control  power  even 
at  reduced  and  zero  rotor  speed.  Not  only  is  this  required  in 
an  actual  aircraft,  the  partial  derivative  matrix  in  the  trim 
section  will  be  singular  if  this  condition  is  not  met  by  the 
control  system  linkage  specified  by  the  input  data.  This 
usually  means  that  the  "fixed  wing”  control  surfaces  must  be 
active  before  the  rotor  is  stopped  and  folded. 

Stopping  the  TFTA  rotor  is  done  most  easily  by  reducing  the 
power  supplied  (J  =  13  card)  and  by  increasing  the  collective 
pitch  of  the  rotors  (J  =  1  card)  to  the  feathered  position. 

The  high  pitch  setting  required  to  feather  (i.e.,  zero 
torque  at  zero  RPM)  the  rotor  can  be  approximated  by  computing 
the  collective  pitch  required  to  set  the  angle  of  attack  at 
the  0.75  blade  radius  equal  to  zero.  Depending  upon  the  rotor 
hub  design,  the  rotor  speeds  at  which  the  flapping  stops 
(J  =  28  card)  are  activated  and/or  the  hub  spring  rate  CJ  = 

25  card)  is  changed,  must  be  specified.  If  the  rotor  must 
be  stopped  at  a  specified  azimuth  location,  a  rotor  brake 
must  be  simulated  CJ  =  18  card).  Once  stopped,  the  TFTA 
rotor  blades  are  folded  aft  uniformly  CJ  =  6  card).  Any  con¬ 
trol  motions  required  during  this  maneuver  must  also  be  speci¬ 
fied  by  appropriate  "J"  cards. 

Stopping  the  HSF  rotor  requires  use  of  the  "J"  cards  mentioned 
above  except  that  the  J  =  27  card  is  used  to  fold  the  blades 
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instead  of  the  J  =  6  card.  Prior  to  stopping  the  rotor, 
the  trim  flight  condition  must  be  determined  which  results 
in  rotor  thrust  and  flapping  near  zero.  This  condition 
minimizes  loads  and  angular  motions  of  the  aircraft  during  the 
stopping  and  folding  maneuver.  Once  the  rotor  is  stopped, 
the  blade  azimuth  position  of  each  blade  must  be  specified 
as  a  function  of  time  by  use  of  as  many  J  =  27  cards  as 
required  to  define  the  relationship.  Blade  feathering  posi¬ 
tions  during  the  fold  cycle  continue  to  be  related  to  the 
rotor  control  positions  by  means  of  the  sine  and  cosine 
functions  of  the  azimuth  angle  as  well  as  the  collective 
pitch.  Minimizing  rolling  and  pitching  moments  during  the 
fold  cycle  currently  requires  several  trials  with  high  fuse¬ 
lage  inertias  which  use  different  rotor  control  relationships 
for  each  trial. 


SECTION  VI 
CONCLUSIONS 

The  analysis  techniques  described  herein  provide  a  versatile  and 
powerful  tool  which  can  provide  data  concerning  performance, 
stability  and  control,  blade  airloads,  and  maneuvering  flight  in 
one  integrated  computer  program. 

The  addition  of  equations  representing  the  stop-fold  rotor  air¬ 
craft  in  the  mathematical  model  enables  accurate  computation  of 
the  characteristics  of  this  concept.  Both  the  horizontal -stop- 
fold  rotor  and  the  tilt-f orward-trail-af t  rotor  characteristics 
may  be  determined. 

The  generality  of  the  mathematical  model  of  the  computer  pro¬ 
gram  enables  the  user  to  simulate  fixed-wing  as  well  as  rotary¬ 
wing  aircraft.  The  former  may  be  jet  or  propeller  driven  and 
the  latter  may  utilize  single,  tandem,  or  side-by-side  rotor 
configurations.  The  generality  of  the  program  requires  over 
400  values  to  be  assigned  to  input  data  to  completely  describe 
the  aircraft  being  simulated. 

The  maneuver  simulation  enables  the  effect  of  changes  of  design 
parameters  to  be  determined.  Flying  qualities  specifications 
which  are  written  assuming  determination  of  compliance  by  flight 
test  may  be  examined  by  simulating  the  required  maneuver  with 
the  computer  program. 

Aerodynamic  interference  effects  are  significant  to  the  accurate 
evaluation  of  stability  and  control  characteristics  and  have 
been  included  in  the  computer  program.  Nonlinear  aerodynamic 
characteristics  can  adversely  affect  the  ability  of  the  computer 
program  to  compute  the  trim  flight  condition. 

The  consistent  use  of  vector  algebra  and  the  complete  matrix 
form  of  angular  rotations  eliminates  all  small  angle  assump¬ 
tions  in  the  analysis. 

It  is  recommended  that  this  computer  program  be  developed 
further  to : 

-  increase  the  sophistication  of  the  mathematical  model 

-  provide  additional  output  data  in  more  convenient  form 

-  include  aeroelastic  blade  loads  if  computer  size  permits 

-  include  additional  equations  describing  automatic  stabi¬ 
lization  functions 

-  compute  linear,  coupled  stability  characteristics 

-  represent  control-free  stability  characteristics  directly 
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Figure  1.  -  Major  Computational  Functions  Flow  Chart 
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Figure  2 .  The  Effect  of  Increment  Size  on 
Stability  Derivatives 
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7  94.1*6 

0.0 

9.0 

0.9 

0.0 

0.0 

FLARRING  ACCEL  MR 7  MAST 

■*74.44? 

-■■74.44? 

0.0 

0.0 

0.9 

9.0 

0.0 

FLARRING  VELOCITY  «RT  HAST 

*.41) 

-*.41? 

0.0 

0.0 

C.O 

0.0 

0.0 

flARRING  LOCATION  MAT  MAST 

1.777 

1.77* 

0.0 

0.0 

C.O 

0.0 

0.0 
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TABLE  V 

AIRLOADS  SUBROUTINE  OUTPUT  DATA 


aiimuth  o- shaft  v-shaft  u-shaft 

MX  S THIBET A1  COSMiTAI  SETA  DOT 

m. 

o 

o 

■M 

o 

28595  -14 

.04529  -25. 

o 

8T 

8* 

m 

*1 

95410  0. 

01491  0. 

99919  -16. 

7(941 

AAO.STA. 

FHI 

ALFHA 

Cl 

CO 

MACH 

LOCAL  VI  IOC.  LAHRoA 

UT 

UF 

1.00000 

0.29S17 

5.41272 

0.50055 

0.01521 

0.4(814 

4. 11 101 

-14.16117 

545.44118 

2.81008 

0.9S000 

0.15811 

5.77568 

0.52592 

0.01568 

0.45491 

6.11101 

-11.91126 

508.15721 

1.40246 

0.90000 

-0.00061 

6.11692 

0.54911 

0.01617 

0.42151 

6.31101 

-11.49915 

470.(5107 

-0.00515 

O.ISOOO 

-0.1(471 

6.41084 

0.57015 

0.01667 

0.1(811 

6.1II0I 

-11.06701 

411.  54460 

-1.41277 

0.10000 

-0.407(2 

6.70972 

0.58851 

0.01714 

0.15474 

6.11101 

-12.61491 

196.21851 

-2.82041 

0. TSOOO 

-0.6  7480 

6.94266 

0.60296 

0.01755 

0.12114 

6.11101 

-12.20282 

158.91217 

-4.22801 

0.70000 

-1.001(6 

7.11169 

0.61245 

0.017(5 

0.2(798 

4.11101 

-11.77072 

121.62622 

-5.61567 

0.45000 

-1.41907 

7.19848 

0.61505 

0.01796 

0.25462 

6.11101 

-11.11861 

2(4.11982 

-7.04111 

0.40000 

-I.9S944 

7.15809 

0.40760 

0.01 779 

0.22127 

4.11101 

-10.906  49 

247.01171 

-8.45092 

o.ssooo 

•2.69IS4 

4.92601 

0.58464 

0.01717 

0.18795 

4.11101 

-10.47418 

209.70758 

-9.(5854 

o.soooo 

-1.71888 

6.17867 

0.51597 

0.01590 

0. 15447 

6.11101 

-10.04228 

172.40118 

-11.26618 

0.45000 

-S. 35946 

5.25(09 

0.44019 

0.01168 

0.12148 

4.11101 

-29.61017 

115.09501 

-12.671(0 

0.40000 

-8.19415 

2*92140 

0.24407 

0.01044 

0.0(845 

4.11101 

-29.17805 

97.7(888 

-14.08142 

O.ISOOO 

-14.16418 

-2  .  74664 

-0.22889 

0.01027 

0.05589 

4.11101 

-28.74594 

60.482  71 

-15.4(901 

0.10000 

•  16.09161 

-21.97606 

-0.85969 

0.10107 

0.02540 

6.11101 

-28.11184 

21.17651 

-16.(9665 

0.29000 

-127.46579 

-115.04821 

0.87075 

1.72244 

0. 02070 

4.11101 

-27.(8171 

-14.12970 

-18.10429 

0.20000 

-159.01154 

-145.91197 

1.01206 

0.65751 

0.04911 

4.11101 

-27.44962 

-51.41591 

-19.71190 

O.ISOOO 

-166.61111 

-152.995  71 

0.92028 

0.41090 

0.0(167 

4.31101 

-27.01752 

-88.  74214 

-21.11954 

0.10000 

-169.(6711 

-155.74956 

0.86547 

0.15217 

0.11441 

4.11101 

-26.58540 

-126.04817 

-22.52718 

o.osooo 

-171.44428 

0.0 

0.0 

0.0 

0.147(1 

4.11101 

-26.15129 

-161.15458 

-21.91478 

A21NUTH 

U- SHAFT  V- 

SHAFT  H-SHAFT 

XX  S1NISE7AI  C0SKE7AI  BETA 

DOT 

91.62015  200. 

28595  -16 

.06529  -25. 

72(55  0. 

95410  0. 

o 

i r 

o 

.99(19  14.78944 

AAO.STA. 

FHI 

ALFHA 

CL 

CO 

HACh 

LOCAL  VI  LOC.  LAM80A 

ut 

OF 

1.00000 

-4.10591 

0.52009 

0.06688 

0.04169 

0.(4414 

6.31101 

•11.90814 

946.01911 

-71.21215 

0.95000 

-4.14594 

0.985  11 

0.12260 

0.02412 

0.81542 

4.11101 

-11.49(99 

908.77026 

-68.98418 

0.90000 

-4.17891 

1.44714 

0.17114 

0.014(0 

0.7(250 

6.31101 

-11.08961 

871.50  1  46 

-66.71599 

0.85000 

-4.42029 

1.90578 

0.21581 

0.01419 

0.74908 

4.31101 

-12.68027 

814.21242 

-64.48782 

0.80000 

-4.46551 

2.16954 

0.25742 

0.011(2 

0.71564 

6.11101 

-12.27090 

796.96118 

-62.21965 

0.75000 

-4.51515 

2.81041 

0.24656 

0.01141 

0.48224 

4.31101 

-11.86156 

759.69414 

-59.49149 

0.70000 

-4.54992 

1.256  15 

0.11164 

0.01140 

0.64(82 

6.11101 

•11.45219 

722.42524 

-57.74112 

0.45000 

-4.41064 

1.69541 

0.14892 

0.01170 

0.41540 

6.11101 

-11.04281 

685.15625 

-55.49511 

0.40000 

•4.49(12 

4.12775 

0.40261 

0.01191 

0.5(148 

4.11101 

-10.41148 

447.(8721 

-51.24496 

0.55000 

•4.77424 

4.551(1 

0.41480 

0.01422 

0.54(56 

4.11101 

-10.22412 

410.61(14 

-50.  44874 

0.50000 

•4.(4004 

4.94401 

0.44554 

0.01462 

0.51515 

4.11101 

-29.(1474 

571. 14412 

-4(. 75061 

0.45000 

-4.95771 

5. 14814 

0.444(7 

0.01 S0( 

0.4(173 

4.31101 

•24.40541 

516.08012 

-46.50244 

0.40000 

-5.04997 

5.75404 

0.52241 

0.01540 

0.44(12 

6.11101 

-28.49605 

49S.(||2( 

-44.25427 

0.15000 

-5.20010 

4.12574 

0.54(68 

0.01614 

0.41441 

6.31101 

•28.58469 

461.54224 

-42.00610 

0. 10000 

-5.15145 

4.47261 

0.57271 

0.01471 

0.1(150 

6.11101 

-28.17714 

424.27144 

-14.75741 

0.25000 

•5.51600 

4.79004 

0.54411 

0.01710 

0.14(04 

6.11101 

-27.76797 

1(7.00419 

-17.50477 

0.20000 

-5.75711 

7.04(75 

0.41278 

0.017(1 

0.11444 

6.11101 

•27.15841 

144.71515 

-15.26158 

0.15000 

•4.01117 

7.294(4 

0.42704 

0.01(27 

0.2(124 

4.11101 

-24.94924 

112.44611 

-11.01141 

0.10000 

-4. 17 ((0 

7.44724 

0.41542 

0.01(55 

0.24791 

6.11101 

•24.51940 

275.14727 

•10.74524 

0.05000 

•4.(1441 

0.0 

0.0 

0.0 

0.21451 

4. 11 1C  1 

•24.11054 

217.92(19 

-2(.  51704 
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TABLE  VI 

VECTOR  ANALYSIS  SUBROUTINE  OUTPUT  DATA 


BEIL  HELICOPTER  |M  360/  MOCHA*  ASA  JO  1 
HELICOPTER  MGIO  BODY  DYNAMICS  ANALYSIS 
COMPILED  10/30/64 
COMPUTEO  0  3/06/70 


HUEY 

CUT  PUT  FIGURES  POA  PINAL  REPORT 
Y«  140  KTS  GW  4S00 


CG-  200 


LEAS!  SGURALS  CUAVE  Pll  START  ING  APTEH  3.644  SLCONOS  MANEUVER  TIME 
Mil  •  AMPL I  T10E*SIAICMEGA*T  4  PMASL  ANCLE)  «  CONSTANT  MITH  OMEGA  ■  C.5C0  CPS 


VAK1AUIE 

AMPLITUOE 

PHASE  ANGLE  IDFGREESI 

CONSTANT 

COEF  OF  COB A 

0  VELOCITY,  IPPl ,  UEG/SEC 

I2.S44 

14B.3S 

0.47288 

0.  66648 

P  VELOCITY,  TPPi  ,  OEC/SEC 

3. 63  34 

4.1134 

-0.14443 

0.60467 

ACT CR1  IHRLSI,  LB 

4331.8 

133.01 

-B4.087 

0. 44826 

P/A  FLAPPING,  NASTI/TPPI,  CFG 

2.8144 

174.13 

-D.  446  07 

0.44417 

LATERAL  FLAPPING,  PASI1/TPPI,  OEG 

0.24311 

139.62 

0.B4244E-01 

0.42461 

0  VELOCITY,  BUOY  AXES,  OEG/SEC 

IT  .442 

131.02 

0 .46444 

0.44434 

BCOY  PITCH  MKT.  (LIGHT  PATH,  OEG 

4.I20R 

70.144 

-11.621 

0.44471 

0  VELUL  1  TV i  1PP I,  OEG/SEC 
RU1UM1  T  HR  LSI  i  LI) 

F/A  FLAPPING,  MASII/TPPI,  OEG 


AMPLITUDE  AND  PHASE  ANGLE  COMPARISONS 

VARIABLES 

/  BOON  PITCH  HAT.  PLIGHT  PATH,  OEG 
/  EOOV  PITCH  NAT  •  PLIGHT  PATH,  DEG 
/  BCOY  PITCH  NAT,  FLIGHT  PATH,  OEG 


AMPLITUDE  RATIO  PHASE  ANGLE  OtFFEAENCF 


LATERAL  FLAPPING,  PASTI/TPP1,  OEG  /  BCOY  PITCH  MAT.  PLIGHT  PATH,  OEG 


3.IIRO 
1044  .7 
0.61323 
0.S4681E-01 


BB. 206 
62.063 
10B.  44 
64.471 


VARIABLE  'A'  AS.A  LINEAR  COMBINATION  OP  VARIABLES  •  B*  ANO  'C*. 

A  «  KBPS  »  KC*C  *  RD 

E  NAME  COEFFICIENT 


P/A  FLAPPING,  MASII/TPPI,  OEG 
0  VELOCITY,  BUOY  ARES,  CEG/SEC 
BODY  PITCH  MKT.  PLIGHT  PATH,  OEG 

CONSTANT 


0.16434 

•O.SB232 

•7,4424 
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TABLE  VII 


FUSELAGE  EQUATIONS 


Aerodynamics 

Lf  =  q*  (XFS20  +  XFS2I  af) 

Df  =  q'  (XFS22  +  XFS23  af  +  XFS24  af2 

+  XFS25  pf2  +  XFS32  cos3pm) 

Yf  =  q'  |  XFS26  +  XFS27  pf  +  XFS28  pf2) 

Mf  =  q*  [(XFSIS  +  XFS16  af ) 

+  XFS32  cos3 p  (XFS31-XMR10)/12 
rm  1 

Nf  =  q'  ( XFS17  +  XFS18  pf ) 

In  addition  to  the  above  fuselage  moments,  the  fuselage  forces 
will  contribute  to  moments  if  the  fuselage  aerodynamic  center 
is  different  from  the  center  of  gravity  by  acting  through  the 
appropriate  moment  arms.  The  following  equation  indicates  this 
effect  on  roll  moment  of  the  fuselage  side  force,  Yf. 

Lf  =  Yf  (XFS04  -  XFS07)/12 

Center  of  Gravity** 

BL  =  XFS06 

SL  =  XFS05  +  3§?§§f-[  (XMR10  -  XFS31)  sinpm 

+  CXMR08  -  XFS30)  (  1  -  cospj] 

WL  =  XFS07  +  [  CXMR10  -  XFS31)  (l  -  cospm) 


-  CXMROG  -  XFS30)  sinpm 


NOTE:  See  Volume  II  for  explanation  of  symbols. 

**Assumes  rotor  masts  tilt  together  and  the  c.g.  buttline 
does  not  vary  with  mast  tilt. 
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TABLE  VIII 


AIRFOIL  CHARACTERISTICS  INPUT  DATA 


Word  Parameter  Units 


Card  A  YXX*  I  Critical  Mach  Number  for  drag 

divergence  for  a=0  Nondim 

2  Mach  number  (Ms)  for  lower  boundary 

of  supersonic  region  Nondim 

3  Maximum  CL,  normal  flow,  M=0  Nondim 

4  Coefficient  of  Mach  number  (M)  in  Nondim 

ClMAX  Equation,  normal  flow 


5  Coefficient  of  M2  in  Equation, 


Normal  Flow  Nondim 

6  Coefficient  of  in  Cjw.y  Equation, 

normal  flow  Nomdim 

7  Maximum  Cj^  reverse  flow,  M=0  Nondim 

Card  B  8  Coefficient  of  Mach  number  (M)  in 

CtMAX  Equation,  Reverse  Flow  Nondim 

9  Coefficient  of  M2in  C^.  Equation 

Reverse  Flow  1«AX  Nondim 

10  Coefficient  of  M^  in  ClwaX  Equation, 

Reverse  flow  Nondim 

11  Tail  boom  bending  coefficient  for  Rad/lb 

correction  of  CL  of  a  fin  or  of  Lifl 

elevator 

12  Drag  coefficient  for  q=0  and  M=0  Nondim 


13  Coefficient  of  a  in  nondivergent  /deg 

drag  equation 

14  Coefficient  of  a2  in  nondivergent  /deg2 

drag  equation 

lard  C  15  Coefficient  used  in  supersonic  drag  Nondim 

equation 
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TABLE  VIII  -  Continued 


Card  C  16 

Maximum  drag  coefficient  for  M<1;  also 
value  of  drag  coefficient  for  1<M<MS 

Nondim 

17 

2  dimensional  lift  curve  slope  for  M=0 

/deg 

18 

Aspect  Ratio 

Nondim 

19 

Not  Used 

- 

20 

Not  Used 

- 

21 

Not  Used 

- 

*  Note:  YXX  represents  the  word  for  the  appropriate  surface; 

ie :  such  as  YFN  for  fin,  YWG  for  wing,  YMR  for  main 
rotor. 
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TABLE  X 


EQUATIONS  OF  MOTION 


Flight  Path 


=  m(U  +  qW  -  rV) 

=  m(V  +  rU  -  pW) 

=  m(W  +  pV  -  qU) 

=  V  -  txz*  *  ’r(Iz  -  V  -  ^ 

=  Iyq  ♦  rp(Ix  -  I2)  +  Ixz(p2  -  r2) 

=  M  -  K,P  *  N(I„  -  I  )  +  ivzqr 


=  q  cos<t>  -  r  sin<t> 

=  p  +  q  sin<t>  tan0  +  r  cos<J>  tan© 
=  (q  sin<f>  ♦  r  cos4>)/cosP 
pcosijf  cos0  cos\|/  sin©  sin<t> 


_  Isinij'  cos0 


-  sin© 


-  sin\J/  cos4> 

sin\J/  sin0  sin<|> 
+  cos\|/  cos4> 

cos0  sin<t> 


cos\| i  sin0  cos<t> 
♦  sinty  sin<|) 

sinij/  sin0  cos<t> 
-  cos\|/  sin<|> 

cos0  cos<J> 


Rotors 


Torque  Supplied  =  Itorsional®  ♦  Torque  Required 
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